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SUMMARY 


The  work  summarized  in  this  report  covers  the  sixth  quarter  of  a  program  with  a 
goal  that  is  threefold:  first,  to  study  the  properties  of  native  point  defect  in  infrared 
focal-plan  array  (IRFPA)  active  and  substrate  materials;  second,  to  study  the  prop¬ 
erties  of  native  point  defects  in  two  classes  of  photonic  materials,  the  wide-gap  II- VI 
compounds  (ZnSe  as  the  prototype  for  which  impurity  properties  will  also  be  calcu¬ 
lated)  and  the  nonlinear  optical  materials  (LiNbOa  as  the  prototype);  and  third  to 
study  the  properties  of  HgZnTe  as  a  VLWIR  detector  material.  Our  accomplishments 
in  the  sixth  quarter  include 

•  Completion  of  an  extensive  paper  on  native  point  defects  in  Hgo.8Cdo.3Te 
as  a  function  of  processing  conditions  for  submission  to  Physical  Review 
B. 

•  Calculation  of  the  binding  energy  of  a  mercury  vacancy  tellurixim  anti¬ 
site  defect  complex  in  Hgo.8Cdo.2Te  and  estimates  of  the  complex  den¬ 
sity  and  its  consequences  on  materials  processing  and  device  perfor¬ 
mance. 

•  Completion  of  the  defect  formation  energies  m  CdTe,  including  gradient 
corrections  to  the  local  density  approximation. 

e  Completion  of  the  defect  formation  energies  in  ZnSe,  including  gradient 
corrections  to  the  local  density  approximation,  and  estimates  of  the 
neutral  defect  concentrations. 

•  Continued  calculation  of  the  defect  ionization  energies  m  CdTe,  ZnSe, 
and  LiNbOa. 

•  Preliminary  prediction  of  the  defect  densities  in  z  =s  0.17  LWIR  Hgi_xZnxTe. 

•  Continued  development  of  a  method  to  calciilate  the  temperature  de¬ 
pendence  of  the  semiconductor  bandgaps,  with  initial  agreement  at 
300  K  both  for  Hgo.78Cdo.33Te  and  GaAs. 

•  Modification  of  our  thermodynamic  codes  for  the  LiNbOs  problem  so 
that  the  stoichiometry  and  temperature  may  be  specified  and  the  defect 
densities  determined. 
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1.  NATIVE  POINT  DEFECTS  IN  HgCdT*  AND  RELATED  IR  MATERIALS 

1.1  NATIVE  POINT  DEFECTS  IN  HfoaCdo  aTt 

We  have  completed  an  extensive  paper  on  the  native  point  defects  in  HgojCdo.aTe; 
a  copy  of  this  paper  is  included  as  Appendix  A. 

1.2  MERCURY  VACANCY  TELLURIUM  ANTISITE  COMPLEX 

We  have  completed  a  preliminary  calculation  of  the  mercury  vacancy  tellurium 
antisite  electronic  binding  energy  and  find  an  energy  of  1.1  eV  per  pair.  This  binding 
energy  corresponds  to  the  reaction 

VkfTe  +  TeHjTe  — >  (VkjTe  +  TeHjTe)patr.  (1) 

From  the  density  of  the  isolated  mercury  vacancy,  Vkg,  and  tellurium  antisites,  XeHg, 
the  density  of  the  pair,  (V^gTe  +  TeHgTe)pa«>,  can  be  calculated  via  the  mass  action 
relation 

[  VngTe  +  T^gTe]pa,> 

If  we  assume  that  the  defect  formation  free  energy  from  changes  to  the  vibrational 
spectrum  are  similar  for  the  isolated  defects  and  the  pair,  we  can  calculate  the  reaction 
constant  as 

=  12H-»exp  (i^)  (3) 

where  [Is]  is  the  density  of  lattice  sites.  The  factor  of  12  is  the  number  of  second 
nearest  neighbor  sites  (but  first  nearest  cation  neighbor  sites)  and  accotmts  for  the 
degeneracy  in  the  number  of  ways  a  pair  can  be  formed  in  which  one  of  the  pairs 
occupies  a  particular  lattice. 

Applying  the  reaction  constant,  we  have  calculated  the  number  of  pairs  present 
in  material  equilibrated  at  500°C  for  the  mercury-rich  and  tellurium-rich  conditions. 
Results  are  summarized  in  Table  1.  We  have  also  estimated  the  concentration  of 
tellurium  antisites  present  following  a  low-temperature  (250*0)  mercury-saturated 
aimeal,  assuming  that  the  mercury  vacancies  can  reach  equilibrium,  but  that  the 
tellurium  antisites  are  frozen-in  at  their  500*0  levels.  The  freezing-in  of  the  anti¬ 
sites  is  proposed  because  of  the  slow  diffusion  rate  expected  for  this  defect  at  lower 
temperatures;  in  this  example  we  have  assumed  that  the  tellurium  antisite  diffusion 
effectively  stops  at  500*0,  and  therefore  the  antisites  are  frozen  in  at  their  500*0 
levels.  In  contrast,  if  the  tellurium  antisites  are  assumed  to  fully  equilibrate  at  250*0 
their  concentrations  would  only  be  ~  6  x  10®  cm“®.  At  this  point  we  do  not  know 
at  what  temperature  the  tellurium  antisite  diffusion  effectively  stops;  this  will  be  de¬ 
termined  by  the  concentration  of  the  vacancy-antisite  pairs  if  the  antisite  diffusion 
is  via  a  vacancy  mechanism.  An  estimate  of  the  critical  temperature  at  which  the 
tellurium  antisite  diffusion  ceases  and  the  frozen-in  tellurium  antisite  concentrations. 


—  Kpair' 


[Vk,Te][TeH,Te] 
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[TcHg],  will  be  made  next  quarter.  If  [Teng]  ^  10^^ - 10^®  cm~®,  this  will  be  an  in¬ 

dication  that  the  antisite  is  the  residual  donor  and  will  allow  us  to  develop  strategies 
for  reducing  the  residual  donor  density. 


Table  1.  Neutral  defect  concentrations  including  the  mercury  vacancy  tellurium 
antisite  pairs.  For  SOO^C,  full  equilibrium  is  assumed  achieved.  For  the  250"C 
mercury-saturated  anneal,  the  total  tellurium  antisite  concentrations  are  assumed 
frozen-in  at  their  500*C  levels. 


Annealing  condition 

iQsniQiiBa 

Emssmsssm 

T=500®C;  Te-saturated 

9.62  xl0“ 

2.58  xlO^* 

T=500®C;  Hg-saturated 

1.37x10^* 

7.52  xlO^® 

T=250*C;  Hg-saturated  following 
a  T=500*C;  Te-saturated  anneal 

3.58x10® 

3.00  xlO^^ 

T=250°C;  Hg-saturated  following 
a  TsSOO^C;  Hg-saturated  aimeal 

3.58x10® 

1.37  xlO” 

1.3  NATIVE  POINT  DEFECTS  IN  CdTs 

We  have  completed  the  calculation  of  the  electronic  contribution  to  the  native 
point  defect  formation  energies  in  CdTe  using  the  full-potential  local-density  ap¬ 
proximation  with  the  linearized  mufhn-tin  orbital  (LMTO)  basis.  Gradient  correc¬ 
tions  were  also  computed.  These  two  contributions  defect  formation  energies  are 
listed  in  Table  2.  These  energies  can  be  used  to  predict  the  neutral  native  defect 
concentrations;  this  calculation  will  be  completed  next  quarter.  The  corresponding 
numbers  for  HgTe  are  listed  in  Table  3.  One  should  not  try  to  infer  the  primary 
defects  from  the  defect  formation  energies  alone,  without  considering  the  vibrational 
entropy  and  vapor  phase  chemical  potential  contribution  to  the  defect  formation  free 
energies. 
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Table  2.  Neutral  native  point  defect  formation  energies  in  CdTe.  LDA  is  the 
contribution  from  the  local-density  calculation  and  GC  is  the  correction  from  the 
gradient  correction  energy. 


Defect 

Defect  reaction 

LDA 

GC 

Total 

Vcd 

CdTe  ^  VcdTe-i-Cd/,.. 

4.152 

-0.463 

3.689 

Vt, 

Cd/,ee  - ►  CdVre 

0.870 

-1-0.315 

1.185 

Cdxe 

2Cdyfee  *  CdCd*!* 

-1.212 

-1-1.091 

-0.121 

Tecd 

2CdTe  — ►  TecdTe-|-2Cd/„, 

7.062 

-1.150 

5.912 

Cd/e. 

Cd^,je  ►  Cd/gj 

0.384 

0.450 

0.834 

^d/j. 

Cdyfee  - *  ^d/j. 

5.424 

-0.698 

4.726 

Te/cd 

CdTe  — ►  Te/cj  -f-  Cdfree 

0.546 

0.391 

0.937 

1.4  TEMPERATURE-DEPENDENT  BANDGAP 

When  the  temperature  of  the  lattice  is  increased,  the  lattice  expands.  This  ex¬ 
pansion,  usually  called  dilation,  decreases  the  bandgap.  The  number  of  phonons  in 
the  lattice  increases  with  temperature.  The  change  in  potentials  caused  by  the  dis¬ 
placement  of  atoms  from  their  equilibrium  positions  introduces  an  electron  phonon 
interaction.  These  interactions,  viewed  as  a  perturbation,  will  change  the  electronic 
states,  that  is,  the  band  structure.  A  major  contribution  to  the  bandgap  change  with 
temperature  arises  from  this  interaction. 

The  electron  acoustic  phonon  interactions  are  usually  parzuneterized  in  terms 
of  wave- vector-  and  energy-independent  drformation  potential  coupling  constants. 
These  constants  are  extracted  from  hydrostatic  experiments  that  measure  the  change 
of  symmetry  point  energies  with  pressure.  In  addition,  only  a  few  phonons,  such  as 
r,  X,  or  L,  are  used;  the  measured  atomic  displacements  due  to  other  phonons  are 
virtually  unknown.  Although  these  calculations  based  on  this  procedure  yield  cor¬ 
rect  trends  in  large-gap  materials,  the  uncertunty  in  measured  vedues,  and  hence  the 
input  parameters,  yields  ambiguous  results  for  small-gap  semiconductors. 

Our  calculations  generalized  the  above  approach  considerably.  We  use  accurate 
tight-binding  band  structures  in  the  calculation  of  the  dynamical  matrix  and  electron- 
phonon  interaction.  In  addition,  we  include  contributions  from  all  six  phonon  modes, 
with  wave  vectors  spanning  the  entire  Brillioun  Zone. 

These  calculations  require  knowledge  of  spatial  dependence  of  atomic  potentials  in 
the  bulk.  We  assume  that  attractive  interatomic  matrix  elements  vary  as  r"*  and  the 
Coulomb  repulsive  interaction  varies  as  The  two  unknowns,  m  and  Vq,  are 

adjusted  so  that  total  energy  calculations  yield  an  accurate  equilibrium  bond  length 
and  bulk  modulus. 

Once  the  spatial  dependence  of  potentials  is  known,  the  dynamical  matrix,  which 
is  the  second  derivative  of  the  potential  divided  by  the  appropriate  atomic  mass,  is 
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calculated  and  diagonalized  to  obtain  the  atomic  displacements  for  a  given  phonon 
wave  vector  and  branch. 

Because  the  atomic  displacements  are  small  compared  to  the  bond  length,  the 
change  in  potentials  is  obtained  by  a  Taylor’s  expansion  and  only  the  first  two  terms 
are  retained.  The  change  in  energy  state  at  a  given  electron  wave  vector  is  obtained 
in  a  second-order,  stationary-state  perturbation  theory.  This  procedure  is  usually 
called  a  rigid  ion  approximation. 

Siunmarizing  the  steps  involved  in  the  calculations;  (a)  Choose  the  energy  state 
whose  temperature  variation  needs  to  be  obtiuned.  Get  the  corresponding  electron 
wave  vector,  k.  (b)  Choose  a  phonon  wave  vector,  q.  (c)  Calculate  the  dynamical 
matrix  and  diagonalize  it  to  get  the  atomic  displacements  for  all  six  modes,  (d)  From  a 
knowledge  of  first  and  second  derivatives  of  the  potentials  in  the  local  basis,  calculate 
the  matrix  elements,  <  k|/fi|k  >  and  |  <  k  —  qj/fslk  >  |^.  (e)  From  the  potential 
well  expressions  derived  from  second-order  pertubation  theory,  the  change  in  energy 
is  obtained.  The  calculation  of  energy  gap  variation  requires  repeating  the  above 
procedure  for  conduction  and  valence  band  edge. 

The  calculations  are  carried  for  the  Hgo.rsCdo.jsTe  alloy. 

The  calculated  change  in  the  gap  at  300  K  is  95  meV,  which  is  in  excellent  agree¬ 
ment  with  the  experimental  value  of  98  meV.  We  are  currently  calculating  the  change 
for  a  larger  temperature  range  in  various  semiconductors.  The  results  will  be  provided 
in  the  next  report. 

2.  NATIVE  POINT  DEFECTS  IN  LWIR  HgZnTe 

We  have  estinaated  the  concentrations  of  native  point  defects  in  z  =  0.17  Hgi-cZn^Te 
based  on  our  results  in  i  =  0.2  Hgi_*Cd*Te.  Both  materials  consist  mostly  of  HgTe, 
the  primary  difference  being  the  difference  in  the  lattice  constants  of  the  material.  For 
HgCdTe  (at  any  i)  the  lattice  constant,  a,  is  6.454  A,  while  for  x  =  0.17  Hgi-^Zn^Te 
a  =  6.378;  we  have  used  an  elastic  continuum  model  to  estimate  the  native  defect 
formation  energy  change  due  to  this  difference  in  lattice  constant.  Results  of  this 
calculation  are  given  in  Table  3. 

We  have  calculated  the  total  (neutral  plus  ionized)  native  defect  concentrations  for 
LWIR  HgCdTe  and  HgZnTe  under  mercury-saturated  conditions  and  as  a  function 
of  temperature;  results  are  given  in  Figure  1.  In  this  calculation  we  have  assumed 
that  the  band  structure  and  its  temperature  dependence  are  the  same  in  HgZnTe  and 
HgCdTe.  While  this  is  not  likely  to  be  rigorously  true,  the  experimental  data  base  for 
HgZnTe  is  far  more  limited  and  thus  we  are  forced  to  make  this  approximation.  Once 
we  have  completed  our  calculation  of  the  temperature  dependence  of  the  bandgap, 
we  will  calculate  the  band  structure  properties  for  HgZnTe  rather  than  borrow  them 
from  the  HgCdTe  experiments. 

From  Figure  1  we  see  that  the  mercury  vacancy,  tellurium  antisite,  and  mercury 
interstitial  densities  in  HgZnTe  are  all  approximately  a  factor  of  two  smaller  than  in 
HgCdTe.  This  difference  is  still  subject  to  refinement  and  may  increase.  However 


Defect  concentration  (cm 


Figure  1.  Defect  concentrations  in  LWIR  HgCdTe 
and  HgZnTe  under  mercury-saturated  conditilons 
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at  present  predicted  differences  in  mercury  vacancy  and  mercury  interstitial  densities 
do  not  account  for  the  experimentally  observed  difference  in  the  mercury  diffusion 
rat*s  between  the  two  materials  and  the  superior  bake  stability  of  HgZnTe.  We  thus 
tentatively  conclude  that  this  difference  in  diffusion  rates  is  not  due  to  differences  in 
the  concentration  of  the  diffusion  species,  but  rather  must  be  due  to  differences  in 
the  migration  energy  of  the  vacancies  and  interstitials  in  the  two  materials.  Because 
the  average  lattice  spacing  in  HgZnTe  is  smaller  than  that  in  HgCdTe,  a  smaller 
migration  energy  in  HgCdTe  is  expected. 


Table  3.  Electronic  contribution  to  the  native  defect  formation  energy  in  HgCdTe 
and  Hgo.17Zno.83Te. 


Defect 

HgCdTe 

HgZnTe 

Vh. 

2.09 

2.14 

Vt. 

1.47 

1.52 

TeHg 

0.46 

0.53 

Hgx. 

3.12 

3.17 

Hg/H. 

1.21 

1.27 

®g/T« 

1.12 

1.16 

Tc/k, 

3.38 

3.42 

Te/cj, 

3.37 

3.45 

3.  WIDE-GAP  1 1- VI  COMPOUNDS  (ZnSe  AS  PROTOTYPE) 

We  have  completed  the  calculation  of  the  native  point  defect  formation  energies  in 
ZnSe  using  the  full-potential  local-density  approximation  with  the  linearized  muffin- 
tin  orbital  (LMTO)  basis.  Gradient  corrections  were  also  computed.  Preliminary 
results  using  a  small  basis  set  size  were  given  in  Quarterly  Technical  Report  3.  The 
calculations  completed  this  quarter  use  a  more  converged  basis  set.  The  energies  are 
summarized  in  Table  4.  The  defect  formation  energy  and  energy  coming  from  the 
modifications  in  the  vibrational  modes  of  the  system  were  reported  in  Quarterly  Tech¬ 
nical  Report  3.  Direct  interpretation  of  the  numbers  appearing  in  Table  4  is  difficult 
because  these  energies  do  not  include  the  chemical  potential  of  the  free  Zn  atom, 
which  can  be  substantial;  this  chemical  potential  is  included  in  our  quasichemical 
analysis  of  the  defect  concentrations. 
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Table  4.  Neutral  native  point  defect  formation  energies  in  ZnSe.  LDA  is  the 
contribution  from  the  local-density  calculation  and  GC  is  the  correction  from  the 
gradient  correction  energy. 


Defect 

Defect  reaction 

LDA 

GC 

Total 

Vzn 

ZnSe 

— > 

VznSe  -f  Zn/,ee 

-0.460 

4.629 

Vs, 

Zn/re« 

— ► 

ZnVs. 

-1-0.322 

1.038 

Zns, 

2Zn/re« 

— ► 

ZnZnse 

-1.373 

-t- 1.002 

-0.371 

Sezn 

2ZnSe 

— ► 

SeznSe  +  2Zn/„e 

9.574 

-1.158 

8.416 

Zn/z. 

Zn/,ee 

— » 

Zn/z. 

4-0.318 

1.421 

Zn/„ 

Zn/,ee 

— ► 

Zn/„ 

-hO.258 

1.616 

Se/z. 

ZnSe 

— ► 

Sc/z«  ^  Znyree 

7.545 

-0.748 

6.797 

Se/z. 

ZnSe 

— ► 

Se/j,  -1-  Zn/ree 

-0.854 

8.244 

In  Figure  2  we  show  a  plot  of  the  neutral  defect  concentrations  as  a  function  of  zinc 
pressure  over  the  single  phase  stability  region  for  material  equilibrated  at  600*C,  The 
zinc  vacancy,  V^a,  and  the  selenium  antisite,  Sezn»  are  the  dominant  defects,  although 
even  at  this  elevated  temperature  their  concentrations  are  not  comparable  to  desired 
dopant  concentrations  and  therefore  might  be  dismissed  as  potential  compensating 
defects.  Because  we  have  not  yet  included  the  ionization  states  of  these  defects, 
this  conclusion  must  be  viewed  as  very  tentative.  We  are  currently  working  on  this 
problem. 

Although  we  have  calculated  the  neutral  native  defect  concentrations  at  600"  C, 
in  fact  we  are  interested  in  what  happens  in  the  material  at  the  low  temperatures 
used  during  molecular  beam  epitaxial  (MBE)  growth.  While  the  calculation  of  the 
equilibrium  defect  densities  at  lower  temperatures  is  straightforward,  a  number  of 
complexities  must  be  considered.  First,  the  effective  zinc  and  selenium  pressure  dur¬ 
ing  the  growth  must  be  determined  so  that  the  chemical  potential  of  these  species  in 
the  solid  can  be  specified.  Second,  there  is  a  question  about  whether  MBE  growth 
constitutes  an  equilibrium,  quasiequilibrium,  or  highly  nonequilibirium  process.  De¬ 
pending  on  which  of  these  most  closely  describes  the  experimental  conditions  used, 
we  can  determine  the  applicability  of  an  equilibrium  or  constrained  quasiequilib¬ 
rium  quasichemical  analysis  to  the  problem.  It  is  most  likely  that  MBE  growth  of 
ZnSe  constitutes  a  quasiequilibrium  situation,  in  which  case  there  are  straightforward 
modifications  to  our  present  formalism  that  will  permit  us  to  calculate  the  defect  con¬ 
centrations  under  quasiequilibrium  conditions.  We  have  a  theory  of  MBE  growth  of 
wide-gap  semiconductors  mostly  cnnompleted.  However,  completing  it  and  applying 
it  to  this  problem  is  not  possible  at  our  current  funding  level.  Therefore,  to  proceed 
for  now,  we  will  make  an  educated  guess  at  what  is  happening  on  the  growth  surface. 
Finally,  one  must  consider  the  mobility  of  the  native  defects  in  the  MBE  material  dur¬ 
ing  growth,  both  on  the  surface  and  in  the  bulk,  to  determine  whether  equilibrium 


or  quasiequilibrium  can  be  reached.  It  is  doubtful  that  defects  such  as  the  selenium 
antisite  will  equilibrate  once  they  are  locked  into  the  bulk,  because  the  diffusion  of 
this  defect  is  likely  to  be  complex,  involving  one  or  more  additional  native  defects. 
We  are  prepared  to  address  these  important  issues  in  our  future  work. 

This  quarter  we  spoke  with  Dr.  M.  Haase  at  3M  Corporation  and  Dr.  Diego 
Olego  at  Philips  Laboratories  about  our  work  in  ZnSe  to  explore  areas  in  which  our 
theory  might  be  of  assistance  to  their  experimental  program.  A  summary  nnnof  this 
correspondence  is  included  as  Appendix  B. 

4.  NONLINEAR  OPTICAL  MATERIALS  (LiNbOa  AS  THE  PROTOTYPE) 

We  have  modified  our  quasichemical  programs  that  predict  the  native  defect  con¬ 
centrations  in  a  binary  or  pseudobinary  material  as  a  function  of  the  temperature  and 
chemical  potential,  to  allow  the  specification  of  the  temperature  and  the  material  sto¬ 
ichiometry  as  the  external  parameters.  This  modification  is  necessary  to  calculate  the 
defect  concentrations  in  LiNbOa,  which,  when  grown  congruently,  is  known  to  have 
a  lithium-to-niobium  ratio  of  0.486:0.514.  By  specifying  this  stoichiometry  in  our 
mass  action  formalism,  we  can  then  predict  the  concentrations  of  the  native  defects. 
The  complexity  comes  in  first  choosing  the  candidate  defects  that  might  account  for 
the  nonstoichiometry  and  secondly  to  calculate  the  formation  energy  of  these  defects. 
Several  models  of  defects  to  account  for  the  lithium-to-niobium  ratio  have  been  pro¬ 
posed  in  the  literature,  which  we  have  discussed  at  some  length  in  previous  reports, 
and  &om  which  we  have  chosen  several  candidate  native  point  defects  for  calculation 
of  the  native  point  defect  densities.  Several  of  the  models  involve  extended  defects 
such  as  stacking  faults  which  will  not  be  included  in  our  analysis. 

We  are  currently  working  on  the  determination  of  the  ionization  states  of  the  native 
defects  in  LiNbOa  so  that  we  may  predict  the  density  of  both  neutral  and  ionized 
defects.  Another  ingredient  necessary  in  the  calculation  of  the  ionized  defect  densities 
is  the  intrinsic  reaction  constant  that  depends  on  the  band  structure  properties  of  the 
material.  Unfortunately,  for  LiNbOa  there  is  less  experimental  information  available 
on  the  band  structure  properties  than  in  the  other  systems  we  are  studying  and 
therefore  we  must  rely  on  our  own  calculations  of  the  electron  and  hole  effective 
density  of  states  for  the  determination  of  the  intrinsic  reaction  constant.  We  have 
begun  to  extract  these  parameters  this  quarter.  Once  these  native  defect  ionization 
energies  and  the  intrinsic  reaction  constant  are  determined,  we  can  then  predict  the 
native  defect  concentrations  as  a  function  of  the  lithium-to-niobium  ratio. 

5.  WORK  PLANNED 

In  the  next  quarter  we  will  continue  to  study  the  mercury  vacancy  tellurium  anti¬ 
site  defect  complex  in  HgCdTe.  Strategies  for  improving  the  calculation  of  the  defect 
densities  in  HgZnTe  will  be  investigated.  We  will  continue  our  study  of  defects  in 
CdTe,  ZnSe,  and  LiNbOa  this  next  quarter,  focusing  on  the  defect  ionization  energies. 
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Once  the  defect  ionization  energies  are  calculated  in  CdTe  and  ZnSe,  the  prediction 
of  the  defect  densities  will  be  straightforward.  We  will  extract  the  band  structure 
properties  of  LiNbOa  for  input  into  the  defect  density  calculation.  We  will  continue 
our  work  on  the  temperature  dependence  of  the  bandgap  in  the  semiconductors. 
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First-principles  calculation  of  native  defect  densities  in  HgojCdo.3Te 


M.  A.  Berding,  M.  van  Schilfgaarde,  and  A.  Sher 
SRI  International,  Menlo  Pari,  Cal(fomia  940i5 
(March  2,  1994) 

We  nee  a  qaaiichemical  feimaliam  to  make  quantitative  piedictiona  of  the  native  point  defect 
denaitiee  in  Hgo.1Cdo.3Te.  The  electronic  contribution  to  the  defect  formation  free  energy  ia  cal¬ 
culated  uaing  the  aelf-conaiatent  firat-prindplea  full-potential  lineariaed  muffin-tin  orbital  method 
and  the  local-denaity  approximation.  A  gradient  correction  ia  added  to  the  LDA  reault  ao  that 
abaolute  reference  to  the  chemical  potential  of  the  mercury  vapor  phaae  can  be  made.  A  Green  *a 
function  approach  baaed  on  a  valence  force  field  plna  a  point  Coulomb  model  ia  uaed  to  calculate  the 
vibrational  contributiona  to  the  defect  fm  energy  (both  energy  and  entropy).  We  find  the  double 
acceptor  mercury  vacancy  ia  the  dominant  defect,  in  agreement  with  previoua  interpretationa  of 
experimenta.  The  tellurium  antiaite  ia  alao  found  to  be  an  important  defect  in  thia  material.  Predic- 
tiona  of  the  low-temperature  hole  concentrationa  are  made  aa  a  function  of  annealing  temperature 
and  compared  urith  available  experimenta.  The  order  of  magnitude  of  our  predictiona  agreea  well 
with  experimental  reaulta,  and  diacrepandea  can  be  attributed  to  contributiona  to  the  free  energy 
that  we  have  neglected  or  to  inaccuradea  in  the  intrinaic  reaction  conatant  uaed.  Snggeatioiu  for 
further  experimental  work  are  made. 
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I.  INTRODUCTION 

The  pseudobinary  aemiconductor  alloy  Hgi.cCdcTe  with  z=0.22  is  currently  the  material  of  choice  for  high- 
performance  detectors  in  the  long-wavelength  infrared  (8-14  ftm).  Unlike  other  II- VI  systems,  both  extriiuic  p- 
and  n-type  doping  can  be  achieved  in  Hgo.sCdo.aTe  although  in  as-grown  material  the  electrical  properties  are 
often  determined  by  native  point  defect  concentrations.  The  dominant  defect  is  believed  to  be  a  double-acceptor 
mercury  vacancy^;  post-growth  low-temperature  mercury-saturated  anneals  are  routine  for  the  reduction  of  the 
mercury  vacancy  concentration.  As  in  other  semiconductors,  it  is  more  difficult  to  establish  the  presence  and  identity 
of  neutral  and  compensating  point  defects,  much  less  to  determine  their  concentrations.  Diffusion  measurements^ 
indicate  the  presence  of  mercury  vacancies  and  mercury  interstitiab,  as  well  as  tellurium  interstitials,  although  no 
unambiguous  determination  of  their  densities  can  be  made  from  these  experiments. 

Although  extended  defects  such  as  dislocations  often  appear  to  be  the  performance  limiter  in  current  state-of- 
the-art  Hg0.sCd0.3Te  devices,^*  a  number  of  mysteries  still  persist  that  may  relate  to  native  point  defects.  (1)  For 
operation  at  40  K,  there  is  a  variation  in  RoAj  and  lifetime  among  pixels  with  no  etch  pits,*  indicative  of  spatial 
nonuniformity  in  the  material  that  is  unrelated  to  dislocations.  (2)  An  as  yet  unidentified  donor  limits  the  minimum 
n-type  carrier  concentrations  obtainable  during  a  mercury-saturated  low-temperature  anneeJ  of  the  material.  While 
the  pressure  and  temperature  dependence  of  this  residu^  donor  does  not  appear  to  correlate  with  the  equilibrium 
dependencies  of  any  native  point  defect,*  the  nearly  universal  presence  of  the  donor  in  liquid  phase  epitaxy  (LPE), 
solid-state  recrystalized  (SSR),  and  molecular  beam  epitaxy  (MBE)  materials  and  its  elusive  nature  do  suggest  that 
a  native  point  defect  is  responsible.  (3)  Undoped  LPE  material  that  has  been  subjected  to  a  low-temperature, 
mercury-saturated  anneal  and  nominally  converted  to  n-type  shows  an  anomalously  low  mobility.  One  interpretation 
is  that  it  is  a  consequence  of  interpenetrating  p-  and  n-type  regions,  with  the  Ugh  effective-mass  holes  lowering 
the  measured  Hall  mobilities.'^  If  this  model  proves  to  be  correct,  it  may  well  be  a  native  point  defect  that  causes 
nonuniform  annealing  of  the  material.  To  overcome  this  low  mobility,  a  donor  impurity  is  added  in  concentrations 
above  that  of  the  unknown  residual  donor.  Thus,  to  lower  the  n-doping  to  desirable  levels,  an  understanding  of  the 
origin  of  the  doping,  and  the  low  mobility  and  a  recipe  for  its  elimination  are  neeeded.  (4)  The  identity  of  the  primary 
Shockley-Reed-Hall  (SRH)  recombination  centers  has  not  been  established;  if  they  can  be  correlated  with  native  point 
defects,  strategies  for  their  elimination  can  be  developed.  (5)  MBE  nuiterial  is  often  n-type  as  grown*  and  may  be 
related  to  a  nonequilibrium  population  of  native  point  defects. 

Unraveling  the  matrix  of  usually  indirect  and  often  contradictory  experimental  data  on  the  native  point  defects  in 
semiconductors  is  a  complicated  task,  especially  when  the  defects  may  be  spatially  varying.  Numerous  theoretical 
efforts  have  utilized  first-principles  methods  to  elucidate  the  properties  of  native  point  defects  in  the  group  IV,  HI-V, 
and  II- VI  semiconductors  (see  for  example  Refs.  9-12).  While  these  studies  have  led  to  much  insight  into  the  properties 
of  the  point  defects,  no  quantitative  predictions  of  the  defect  densities  were  made.  Several  earlier  theoretical  studies 
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have  looked  at  the  properties  of  defects  in  HgCdTe^^^*  although  once  again  no  quantitative  predictions  of  defect 
densities  were  given. 

Our  goal  in  this  paper  is  to  theoretically  identify  the  important  native  defects  in  BgCdTe,  to  calculate  their  densities 
as  a  function  of  growth  and  processing  conditions,  to  substantiate  the  experimentally  deduced  properties  of  the  native 
defects,  and  to  begin  to  unravel  the  remaining  mysteries  in  this  material.  To  express  the  concentration  of  the  native 
point  defects  in  terms  of  their  formation  free  energies,  we  employ  the  quasichemical  formalism.  In  addition  to  the 
electron  and  hole,  we  have  included  eight  native  point  defects  (and  their  ionized  species)  in  the  analysis;  the  mercury 
and  tellurium  vacancies,  the  mercury  and  tellurium  antisites,  and  two  types  of  mercury  and  tellurium  tetrahedral 
interstitials  -  one  surrounded  by  four  cation  near-neighbors  and  one  surrounded  by  four  anion  near-neighbors.  As  we 
will  show,  we  have  attempted  to  incorporate  all  of  the  important  contributions  to  the  free  energy  and  adopt  a  first- 
principles  approach  for  most  of  the  quantities  we  calculate.  The  only  significant  empirical  data  we  employ  are  those 
needed  to  obtain  the  temperature-dependent  intrinsic  reaction  constant.  Our  calculated  native  defect  concentrations 
are  in  quantitative  agreement  with  the  available  experiment  data.  To  our  knowledge,  this  is  the  first  attempt  in 
any  material  to  carry  through  a  comprehensive  theory  of  the  native  point  defect  concentrations  essentially  from  first 
principles.'^ 

A  number  of  features  make  our  study  of  defects  in  Eg0.sCd0.3Te  unique,  and  permit  us  to  calculate  absolute  defect 
concentrations: 

(i)  To  calculate  the  electronic  contribution  to  the  defect  formation  free  energies  we  employ  the  self-consistent  first- 
principles  full-potential  linearized  muffin-tin  orbital  (FP-LMTO)  method'*  and  the  local-density  approximation. 
The  LMTO  method  is  well  suited  for  compounds  containing  d  electrons,  such  as  Hgi-cCd,Te. 

(ii)  Because  the  local-density  approximation  (LDA)  overbinds,  we  have  also  employed  gradient  corrections  to  the 
LDA  of  the  Langreth-Mehl-Hu  type.'*  These  corrections  greatly  improve  the  overbinding  found  in  the  LDA.*°  We 
believe  the  calculated  energies  are  precise  enough  that  we  naay  make  comparison  with  atoms  referenced  to  the  free 
atom,  and  therefore  by  combining  these  energies  with  the  translational  energy  of  the  atoms  in  the  vapor  phase,  are 
able  to  calculate  the  chemical  potential  for  a  monoatomic  mercury  vapor. 

(iii)  A  Green’s  function  fortnalism  within  a  valence  force  model  plus  point-charge  ionic  model  is  used  to  calculate 
the  vibrational  contribution  to  the  defect  formation  free  energy,  both  the  enthalpy  and  entropy. 

(iv)  The  combination  of  the  electronic,  translational,  and  vibrational  free  energies  calculated  in  (i)-(iii)  encompass 
the  primary  contributions  to  the  total  defect  formation  free  energies  when  referenced  to  a  mercury  vapor.  The  total 
defect  formation  free  energy  is  then  incorporated  into  a  quasichemical  formalism,*'  and  predictions  of  absolute  defect 
concentrations  as  a  function  of  the  thermodynamic  vanities  temperature  and  pressure  can  be  made. 

The  remainder  of  the  paper  is  organized  as  follows.  In  Section  II  we  describe  the  quasichemical  formalism  used  to 
calculate  the  neutral  native  defect  concentrations  and  its  extension  for  ionized  defects  and  alloys.  The  calculations  of 
the  defect  formation  free  energies  are  discussed  in  Sections  III  and  IV  for  the  electronic  and  vibrational  contributions, 
respectively.  In  Section  V  we  present  the  results  of  our  calculations  and  a  comparison  with  available  experimental 
results.  We  end  with  a  brief  summary  and  conclusions  of  our  work  in  Section  VI. 


II.  DEFECT  CONCENTRATIONS 
A.  Quasichemical  formalism  in  compounds 

We  begin  by  outlining  the  formalism  employed  to  calculate  the  defect  concentrations  in  a  compound  as  a  function 
of  external  parameters.  In  Section  II C  we  discuss  extensions  of  the  formalism  necessary  to  treat  the  low-z  alloy 
Hgi-,Cd,Te. 

The  defect  reactions  for  the  compound  AC  to  be  considered  in  this  paper  are  listed  in  Table  I.  We  have  chosen 
the  AC  unit  ceil  and  A  in  the  phase  R  as  our  reference  states  for  the  calculation  of  the  reactions’  free  energies.  FYom 
Gibbs’  phase  rule  we  know  that  for  a  system  of  two  components  (A  and  C)  and  two  phases  (one  of  which  is  the  AC 
zincblende  solid)  there  are  two  degrees  of  freedom.  For  this  paper  we  shall  assume  that  temperature  is  one  degree  and 
that  the  chemical  potential  of  an  external  reservoir  of  A  or  C  atoms  is  the  other.  The  chemical  potential  is  chosen  to 
be  consistent  with  the  experimental  situation  to  be  modeled.  We  shall  choose  our  reference  state  to  be  the  mercury 
vapor  and  therefore  choose  to  reference  our  reaction  energies  to  the  AC  solid  and  An,  which  will  be  taken  as  the 
monoatomic  mercury  vapor  (extensions  of  this  analysis  to  the  alloys  are  discussed  in  Section  V). 

Reference  to  other  reservoirs  can  be  obtained  by  considering  the  additional  reaction 

AC^Afl  +  Ci,  (la) 


or  for  molecular  species  as  the  reference 
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AC 


(lb) 


fix  flc 

where  n,  is  the  appropriate  integer,  and  by  taking  linear  combinations  of  these  with  the  reactions  in  Table  1. 

We  are  not  restricted  to  specifying  the  temperature  and  a  chemical  potential  of  one  of  the  species  as  the  two  degrees 
of  freedom,  but  could  choose  instead  another  particular  set,  including,  for  example,  the  stoichiometry  of  the  system.^” 
While  other  choices  are  possible,  usually  one  does  not  know  a  phori  the  material’s  stoichiometry. 

Assuming  that  the  defect  densities  are  small  and  that  they  are  noninteracting,  from  the  law  of  mass  action  we  can 
write  the  reaction  constant  for  each  of  the  neutral  defect  reactions  in  Table  I  as^^’^^ 


where  [X’*]  is  the  density  of  the  neutral  defect  X,  0  is  the  number  of  unit  cells  per  volume  and  converts  [X‘‘]  &om 
site  fraction  to  defects  per  unit  volume,  and  fjrx  is  the  free  energy  for  the  neutral  defect  reaction.  The  free  energy 
for  any  defect  X  can  be  written  as  the  sum 

Fx  =  +  /jt  •"*  +  ^bT  b(G),  (3) 

where  ka  is  Boltzmann’s  constant  and  Fx*  is  the  vibrational,  Fx*‘*  is  the  electronic,  and  Fx***'  is  the  translational 
energy  contribution  to  the  reaction  free  energy,  and  G  accounts  for  the  degeneracies  of  the  reactants.  In  the  quasi¬ 
chemical  approximation  (QCA),  G  =  {,9c9d)I{QaSb)  fot  the  reaction  A  +  B  C  •¥  D,  where  9i  is  the  degeneracy 
of  reactant  t.  Once  Fx*  m  know  for  a  given  reaction,  Kx*  can  be  evaluated  and  the  defect  concentration  can  be 
determined.  The  difficulty,  of  course,  is  in  the  evaluation  of  Fx* ,  which  is  discussed  in  Sections  III  and  IV. 


B.  Ionized  defects  and  the  intrinsic  reaction  constant 


The  above  discussion  applies  to  the  neutral  defect  densities.  In  most  semiconductors  the  native  point  defect  will 
have  one  or  more  localized  levels  in  the  band  gap,  allowing  for  multiple  ionization  states  of  the  defect.  We  thus  need 
to  calculate  the  concentration  of  these  ionized  defects,  in  addition  to  the  neutral  concentrations  discussed  above,  to 
obtain  the  total  defect  populations.  Once  the  energies  of  the  localized  levels  are  determined  (Section  III  E)  their 
populations  can  be  calculated  via 


for  an  acceptor  and 


[X'] 

-iilml 

[X-] 

<  kaT 

[X*] 

=  |il.xp| 
9X* 

fEd-HF-rx**  +  n'* 

[X*] 

K  kaT 
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(5) 


for  a  donor  state  of  the  defect  X.  A  bullet  superscript  indicates  a  positive  charge  and  a  prime  a  negative  charge,  £« 
and  Ed  are  the  acceptor  and  donor  one-electron  ionization  energies  with  respect  to  the  val^  '  ce  and  conduction  band 
(both  defined  as  positive  for  states  in  the  gap),  and  /ip  is  the  Fermi  energy.  Although  the  last  term  in  the  exponential, 
-I-  fx'*’  •tould  rightly  be  there,  and  corresponds  to  the  difference  in  the  vibrational  free  energy  of  the  neutral 
and  ionized  defect,  it  has  never  been  considered  previously  and  for  the  present  we  shall  neglect  it,  too. 

For  multiply  ionized  defects  with  positive  Hubbard  U’s  (reactions  7  and  8  in  Table  I),  the  above  expression  for  the 
number  of  ionized  acceptors  generalizes  to 


(6) 


where  z  is  an  integer  and  E^  is  the  ionization  energy  of  the  ionization  level.  A  similar  generalization  applies  for 
the  donor  levels. 

For  each  ionized  defect  concentration,  we  introduced  one  equation.  In  addition,  though,  we  have  two  new  unknowns, 
the  Fermi  energy  and  either  the  electron  or  hole  concentration.  Thus,  two  additional  equations  are  needed. 

First,  we  have  the  additional  reaction  for  the  generation  of  electron-hole  pairs  across  the  band  gap,  reaction  9  in 
Table  I,  and  the  corresponding  intrinsic  reaction  constant 


A-A 


(7) 


=  [*•][«']  =pn 

where  p  s  [h*]  and  n  =  [e'].  In  general,  Kpn  depends  on  the  shapes  of  the  conduction  and  valence  bands,  the  band 
gap  energy,  the  Fermi  energy  (for  degenerate  statistics),  and  the  temperature  variation  of  these  quantities.  Several 
lirnits  are  often  encountered  in  the  evaluation  of  Kpn-  First,  when  the  conduction  and  valence  bands  are  parabolic, 
although  not  necessarily  isotropic,  E  k  k^,  and  the  reaction  constant  can  be  written  in  terms  of  the  Fermi-Dirac 
integrals  as 


Kpn  =  4  (mAm,)*/a 

(®) 

where  is  the  Fermi-Dirac  function;  £«,  E^,  and  pr  the  conduction  band,  valence  band,  and  Fermi  energies, 
respectively;  and  m«  are  the  hole  and  electron  density-of-atates  effective  masses,  respectively;  and  h  ii  Planck’s 
constant.  In  the  nondegenerate  limit,  this  reduces  to  the  familiar  expression 

Kpn  =  4  (mfcm.)»/»exp  (8) 

which  is  independent  of  the  Fermi  energy.  For  the  general  case,  which  will  apply  even  at  moderate  temperatures  for 
narrow-gap  HgCdTe,  Kpn  depends  on  the  extrinsic  carrier  densities  through  its  dependence  on  pp-  Moreover,  the 
conduction  band  in  this  narrow-gap  material  is  not  well  represented  by  a  parabola,  but  is  rather  more  hyperbolic.’^ 
Because  the  hyperbolic  dispersion  relation  approaches  a  linear  dependence  away  from  the  band  edge,  we  will  for 
the  present  assume  a  linear  dispersion  relationship  of  the  form  E  =  ak.  In  this  case  the  intrinsic  reaction  constant 
becomes 


where  is  Fermi-Dirac  integral  of  order  2. 

The  requirement  of  charge  neutrality  leads  to  a  second  additional  equation; 

E  E  *  w'l + = E  E  *  w) + 

it  it 

where  t  sums  over  the  various  defects  and  x  sums  over  the  various  ionization  states  of  the  defect  Xi- 


(10) 


(11) 


C.  Quasichemical  formalism  in  alloys 

We  wish  to  generalize  the  above  formalism  to  the  ideal  cation  substituted  pseudobinary  alloy  Ai-gBgC.  It  is 
perhaps  easiest  to  demonstrate  the  generalization  with  a  specific  defect  reaction,  for  example  a  neutral  vacancy  on 
the  cation  sublattice.  In  the  compound  AC  the  formation  reaction  is  given  by  reaction  1,  Table  I: 

AC^V^C  +  Ar  (12) 

where  indicates  a  vacancy  on  the  cation  sublattice.  In  the  alloy  a  neutral  vacancy  on  the  a  sublattice  can  form 
via  the  same  reaction.  The  only  difference  comes  in  the  evaluation  of  the  reaction  coiutant.  In  the  compound  the 
vacancy  density  is  given  by 

ra=*«p(^)  (13) 

while  in  the  alloy  it  is  given  by 

[V;>‘]  =  «(l-*)exp^^j  (14) 
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where  u  the  reaction  free  energy  correaponding  to  Eq.  (12)  in  the  alloy.  The  factor  of  (1  -  z)  result*  from  the 

coniigurational  entropy  contribution  to  the  chemical  potential  of  A  on  a  lattice  site  (/«),  ib£71n([A(«]/[/«])  s:  kaTlaix). 
Because  we  have  assumed  that  the  defect  concentrations  are  small,  the  coidigurational  entropy  contribution  to  the 
chemical  potential  of  A  in  the  compound  AC  is  ibsTln([i4j(]/[/s])  a  0.  In  both  the  compound  and  the  alloy  9  is  the 
same. 

In  addition  to  Eq.  (12),  in  the  alloy  the  vacancy  on  the  a  sublattice  can  also  form  via  the  reaction 

BC^V^C  +  Bh.  (18) 

with 

[Cl  =  j  (i«) 

where  Fy^  >•  the  reaction  fne  energy  corresponding  to  Eq.  (15)  in  the  alloy.  Now  the  vacancy  concentrations 
predicted  by  Eqs.  (14)  and  (16)  must  be  equal  and  thus 

(l-.)«.p(-^)=(.)«p(^).  (17) 

It  is  ^>parent  that  this  simply  corresponds  to  the  difference  of  Eqs.  (12)  and  (15) 

Ar  +  BC~^AC-¥Br  (18) 

that  is,  the  exchange  of  an  A  and  B  on  a  lattice  site,  which  is  a  reaction  in  the  alloy  in  equilibrium  with  Ar  and  Br. 
An  andysis  similar  to  the  above  applies  for  the  anion  antisite,  Cq,  which  also  substitutes  on  the  cation  sublattice. 

In  the  paeudobinary  alloy,  the  vacancy  free  energies  Fy»  and  Fy^  appearing  in  Eqs.  (13)  and  (14)  may  be  different. 
The  energy  of  a  vacancy  depends  on  the  local  configuration  of  tlie  sunounding  lattice;  this  changes  in  the  second 
and  more  distant  neighbor  shells  for  the  cation  vacancy  and  the  anion  antisite  in  the  cation  substituted  alloys.  A 
completely  rigorous  approach  would  treat  each  kind  of  vacancy  uniquely;  indeed  in  a  previous  work^*  we  found  a 
configuration  dependence  of  the  vacancy  in  the  Aq.sBq.sC  lattice,  of  several  tenths  of  an  eV,  varying  iq>proximately 
linearly  in  the  number  of  A  atoms  in  the  second  neighW  shell.  Here,  we  have  ignored  this  refinement  and  assumed 
the  A  vacancy  and  C  antisite  surroundings  are  totally  of  species  A.  This  is  justified  to  some  extent  because  we  are 
interested  in  low  z  compositions  of  Hgi.cCdcTe. 

For  defects  on  the  C  sublattice  (e.g.,  the  anion  vacancy  and  the  cation  antisites)  as  well  as  interstitial  atoms 
in  certain  tetrahedrrd  sites,  one  must  more  carefully  consider  the  complications  from  the  alloy,  because  disorder  is 
found  already  for  nearest  neighbors.  For  systems  in  which  these  classes  of  defects  are  important,  the  configuration 
dependence  of  the  surrounding  sublattice  must  be  taken  into  account.  For  low  z  Hgi-rCdcTe  we  find  the  densities 
of  these  defects  to  be  quite  low,  and  thus  the  error  incurred  in  using  the  electronic  energies  calculated  for  the  pure 
AC  compound  will  not  impact  the  major  conclusions  of  our  work. 

Finally,  the  band  gap  is  one  other  important  consideration  when  comparing  the  alloy  to  the  pure  compound.  This 
is  of  particular  importance  for  the  intrinsic  reaction  constant  and  we  employed  an  empirical  fit  to  the  temperature 
dependence  of  the  alloy  band  gap  and  intrinsic  carrier  concentration  that  are  used  to  calculate  Kpn-  This  is  presented 
in  detail  in  Section  III  G. 

III.  CALCULATION  OF  DEFECT  FORMATION  FREE  ENERGY:  ELECTRONIC  CONTRIBUTION 

A.  Full  potential  LMTO  calculations 

Total  energy  calculations  for  the  defect  reactions  were  all  obtained  with  a  full-potential  version**  of  the  LMTO 
method,  in  the  local  density  functional  approximation  of  Barth  and  Hedin.**  This  method  has  been  tested  extensively 
for  most  of  the  elemental  ap  and  d  bonded  solids,  the  II- VI,  III-V  and  column  IV  semiconductors,  and  a  host  of  other 
solids.  Our  results,  which  will  be  reported  elsewhere,*®  show  good  agreement  with  experiments  for  all  systems  studied, 
with  small  and  systematic  errors  in  structural  and  mechanical  properties.  The  most  prominent  error,  particularly  for 
the  present  purposes,  is  the  overbinding  of  the  solid. 


A-6 


In  the  FP-LMTO  method,  the  only  important  approximation  we  make  beyond  the  local-deneity  approximation,  liee 
in  the  treatment  of  the  interstitial  matrix  elements.  The  LMTO  method  employs  an  atom-centered  basis,  represented 
by  Hankel  functions  in  the  interstitial.  For  the  calculations  presented  here,  the  basis  consisted  of  a  ‘triple  kappa”  basis 
22  orbitals  per  atom,  with  energies  —0.01,  -1,  and  -2.3  Ry  for  the  s  and  p  orbitals,  and  -0.01  and  -1  Ry  for  the  d. 
Inside  the  muffin-tin  (MT)  spheres,  wave  functions  are  represented  by  spherical  harmonics  and  numerically  tabulated 
radial  functions.  The  electron  density  and  potential  can  be  similarly  represented,  since  the  density  generated  by  a 
Hamiltonian  is  obtained  by  summing  over  the  eigenvectors.  Outside  the  MT  spheres,  another  treatment  is  necessary. 
Methfessel^*  developed  a  simple,  efficient  way  to  represent  the  density  and  potential  in  the  interstitial  by  extrapolation 
from  the  edges  of  MT  spheres,  where  the  value  is  well  known.  The  electron  density  is  represented  in  the  interstitial 
..  as  a  linear  combination  of  Hankel  functions  that  are  chosen  to  match  the  value  and  slope  of  the  function  at  each  MT 
sphere.  Two  Hankels  per  site  and  Im  are  enough  to  match  the  values  and  slopes  at  aU  MT  spheres.  This  representation 
of  the  density  throughout  the  interstitial  is  approximate,  although  it  becomes  exact  near  any  MT  sphere.  Extensive 
tests  show  that  the  approximation  works  very  well  for  close-packed  systems,  but  the  errors  can  become  significant 
when  the  packing  is  poor.  To  ensure  a  good  fit  to  the  charge  dezuity  and  potential  in  the  interstitial  region  of  the 
sincblende  solids,  we  include  empty  spheres  at  each  tetrahedral  interstitial  site  (rendering  the  sphere  packing  bcc  for 
the  ideal  lattice).  In  addition,  we  added  orbitals  to  the  basis  by  centering  them  on  the  empty  spheres.  Addition  of  2s 
and  2p  orbitals  changed  the  energy  by  approximately  0.1  mRy/atom,  showing  that  the  basis  is  nearly  complete. 

To  assess  the  validity  of  the  interstitid  approximation  for  the  representation  of  the  charge  density  and  interstitial 
matrix  elements,  an  alternative  approach  was  developed, which  is  similar  to  a  procedure  described  by  Jones.’^  When 
calculated  in  this  way  the  total  energies  changed  by  approximately  ImRy/atom,  showing  that  the  approximation  is 
a  good  one. 

Both  the  charge  density  inside  the  spheres  and  the  tails  of  Hankel  functions  centered  on  a  neighboring  sphere  were 
expanded  to  1  =  6.  We  estimate  that  the  error  introduced  by  truncation  at  /  =  6  to  be  about  ImRy/atom,  in  line 
with  other  errors  in  the  method.  The  core  was  allowed  to  relax  during  the  self-consistency  cycle.  The  semicore  d 
electrons  in  the  tellurium  were  treated  explicitly  as  valence  states  in  a  second  panel;  explicit  treatment  of  these  states 
was  found  to  introduce  a  small  but  significant  correction  to  the  total  energy.  For  the  16-atom  cells,  the  Brillouin 
zone  integrals  were  done  by  a  sampling  method  for  the  charge  density,  and  the  linear  tetrahedron  method  for  the 
band-structure  energy,  augmented  by  Blbchl  weights,  and  a  mesh  of  four  divisions  was  used  (6  k-points).  Checks 
showed  that  this  was  sufficient  to  converge  the  energy  to  1  mRy/cell. 


B.  Supercell  approximation 

Supercells  are  used  in  which  a  periodic  array  of  defects  is  constructed.  Defect  formation  energies  are  calculated 
from  a  difference  in  total  energies  of  the  compound  with  and  without  the  defect.  For  example,  if  we  denote  CjiV^) 
as  the  energy  of  a  supercell  containing  j  lattice  sites  and  one  A  vacancy,  the  energy  for  defect  reaction  1  in  Table  I  is 
given  by 


=  F;(Vx)  +  E{Ar)  -  SjiAB)  (19) 

where  E{Ar)  is  the  energy  of  an  A  atom  in  the  reference  state  R  and  €j[AB)  =  jE{AB)  where  Eab  i»  the  energy 
of  an  ideal  AB  unit  cell.  For  some  defects  the  number  of  lattice  sites  changes  in  the  reaction;  for  example  for  the 
formation  of  the  B  antisite  via  reaction  4  in  Table  I,  the  formation  energy  is  given  by 


E{Bj,)  =  F;(flx)  +  2E(Ah)  -  £i+2{AB). 


(20) 


Because  we  wish  to  calculate  the  formation  energies  in  the  dilute  limit,  we  use  the  largest  supercell  computationaUy 
feasible.  For  this  paper,  all  calculations  were  done  using  16-atom  supercells. 


C.  Gradient  corrections  to  the  locai  density 

The  local-density  approximation  generally  overbinds  the  solids.  Several  systematic  extensions  of  the  local-density 
function  have  been  proposed  that  are  based  on  generalized  gradient  approximation  for  the  exchange  and  correlation 
energies.  We  have  considered  one  of  these  extensions,  that  proposed  by  Langreth  and  Mehl,^*  and  have  examined 
the  systematica  in  the  gradient  corrections  to  the  lattice  constants,  cohesive  energies,  bulk  modulus,  and  other  elastic 
constants  for  a  wide  array  of  solids^®;  preliminary  results  of  that  work  for  the  zincblende  semiconductors  are  shown  in 
Table  II.  With  few  exceptions,  the  gradient  correction  systematically  improves  the  prediction  of  the  cohesive  energy. 
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although  the  prediction  of  the  elastic  constants  often  worsens  slightly.  The  improvement  in  the  cohesive  energy  is 
largely  due  to  improvement  in  the  calculation  of  the  total  energy  of  the  free  atoms,  rather  than  the  solid. 

Here,  we  are  interested  in  calculating  total  energies  for  reactions  in  which  a  constituent  is  exchanged  between  the 
solid  and  the  vapor,  and  thus  the  errors  inherent  in  the  local-density  calculation  of  the  cohesive  energy  will  be  present 
in  these  energies  idso.  Because  the  relaxations  do  not  change  significantly  when  gradient  corrections  are  added, 
we  have  completed  the  majority  of  the  calculation,  including  the  relaxation,  within  the  FP-LMTO.  The  gradient 
correction  energy,  calculated  at  the  LDA-determined  relaxed  positions,  is  then  added  to  the  LDA  energy. 


D.  Relaxation 

In  general  the  lattice  relaxes  in  the  presence  of  a  defect,  thereby  lowering  the  lattice  energy.  In  the  dilute  defect 
limit,  the  radial  relaxation  of  the  lattice  extends  to  infinity.’*  In  the  supercells  we  account  for  this  relaxation  by 
allowing  the  overall  lattice  constant  of  the  supercell  to  relax  to  mininndze  the  supercell  total  energy.  Because  HgTe 
and  CdTe  are  nearly  lattice  matched,  and  their  elastic  constants  are  the  same,  the  defect  relaxations  in  pure  HgTe 
should  be  comparable  to  those  in  the  HgCdTe  alloys.  Second,  for  the  most  important  defects,  we  permit  the  radial 
relaxation  of  the  defect  near-neighbor  atoms.  For  the  on-site  defects  (the  vacancies  and  antisites)  we  permit  only  the 
nearest-neighbor  atoms  to  relax.  Estimates  of  these  relaxation  energies  are  given  in  Section  V. 

Relaxation  energies  are  calculated  only  for  the  neutral  defects  and  are  assumed  comparable  in  the  ionized  defects. 
Nonradial  relaxations  such  as  the  trigonal  and  tetragonal  distortions  that  split  the  degeneracy  of  the  triply  degenerate 
T]  states  may  be  important  and  may  differ  substantially  for  the  different  charge  states  of  the  system.  Because  the 
symmetry  of  the  distortion  depends  on  the  charge  state  of  the  defect,  distortions  and  charge  states  must  be  treated 
simultaneously.  These  distortions  have  not  been  considered  in  this  paper. 


E.  Localized  defect  levels 

The  calculation  of  the  ionization  states  of  the  defects  is  perhaps  the  most  difficult  part  of  the  calculation  of  the 
native  defect  concentrations,  in  large  part  because  of  the  inadequacies  of  the  LDA  in  predicting  the  band  gap  of 
the  semiconductors.  This  is  additionally  complicated  by  the  fact  th8^  our  calculations  were  done  for  HgTe,  which  is 
known  experimentally  to  be  a  semimetal  with  a  negative  band  gap  of  -0.3  eV,  so  that  even  if  the  LDA  band  gap  were 
correct,  we  would  still  have  a  zero-gap  material.  Furthermore,  because  the  Coulomb  fields  associated  with  a  defect 
may  be  extended,  we  expect  that  very  large  supercells  will  be  needed  to  isolate  the  localized  levels  of  an  individual 
defect. 

We  have  developed  a  method  to  calculate  the  location  of  localized  defect  levels  in  the  band  gap  and  have  applied 
it  to  the  arsenic  antisite  defect  in  GaAs.  This  defect  was  chosen  because  of  its  technological  importance  and  because 
these  levels  have  been  determined  experimentally  by  Weber  et  al.”  Calculations  were  done  within  the  ASA  so  that 
we  could  examine  the  convergence  of  our  results  going  to  large  (128  atom)  supercells.  Our  approach  is  similar  to  that 
discussed  by  van  der  Waal  et  al.*°  in  which  the  shift  in  the  Fermi  level  is  examined  as  electrons  are  added  to  (or 
removed  from)  the  defect,  with  a  compensating  uniform  background  charge  added  so  as  to  maintain  charge  neutrality. 
We  find  good  agreement  with  experiments  of  the  two  antisite  donor  levels.  We  also  have  found  these  energies  agreed 
closely  with  the  positions  of  peaks  in  the  density  of  states,  when  referenced  to  the  top  of  the  valence  band.  Details  of 
the  calculation  will  been  given  elsewhere. 

Because  the  compositions  of  HgCdTe  of  interest  here  have  narrow  band  gaps,  the  determination  of  the  exact  location 
of  the  defect  levels  in  the  band  gap  is  not  as  important  for  the  purpose  of  calculating  the  defect  concentrations  as 
in  a  wider-gap  semiconductor  such  as  GaAs,  although  the  identification  of  the  position  of  defect  levels  is  useful  in 
understanding  mechanisms  limiting  carrier  lifetimes.  We  have  used  the  54-atom  supercells  of  HgTe  to  determine  the 
type  (acceptor  or  donor)  of  the  various  native  defects  based  on  the  position  of  the  Fermi  level  with  respect  to  the 
states  that  lie  within  ~  0.1  eV  above  the  valence  band  edge.  In  addition,  an  assessment  of  whether  the  state  is  a 
single  or  double  donor  or  acceptor  and  whether  the  state  is  shallow  or  deep  has  been  made  based  on  the  position 
of  the  density  of  states  peaks.  For  the  mercury  vacancy,  we  follow  the  arguments  of  Cooper  and  Harrison,*^  and 
assume  that  it  is  a  negative-U  center,  with  the  neutral  and  double  acceptor  states  being  the  only  observable  states; 
this  assumption  is  consistent  with  the  observation  that  the  mercury  vacancies  are  always  found  to  be  doubly  ionized 
acceptors.^ 


F.  Ionization  atate  degeneracy 


The  degeneracy  of  the  various  ionization  states  of  each  defect  may  differ  and  must  be  determined^’  to  complete 
the  calculation  of  the  density  of  ionized  defects.  As  an  example,  we  consider  the  A  vacancy  in  a  II- VI  material,  and 
use  tight-binding  language  for  the  purpose  of  discussing  the  defect  states.  There  are  four  dangling  anion  hybrids, 
each  donating  1.5  electrons  to  the  system,  for  a  total  of  six  electrons  at  the  vacancy  site.  Although  we  have  not 
explicitly  calculated  it,  for  the  purpose  of  computing  state  degeneracy,  we  assume  that  a  symmetry  lowering  Jahn- 
Teller  distortion  will  take  place  whenever  there  is  a  state  degeneracy  beyond  two  (for  spin)  and  a  partial  occupancy 
of  that  state.  Thus,  we  assume  that  the  highest  filled  vacancy  level  in  the  neutral  state  is  doubly  occupied  with  oac 
electron  spin-up  and  one  spin-down,  and  that  the  level  can  accept  no  other  electron.  Because  there  is  only  one  unique 
configuration  for  this  state,  the  state  has  a  degeneracy  of  one.  For  the  single  acceptor  state  in  which  one  electron 
has  been  added  to  the  vacancy,  the  extra  electron  can  either  go  in  spin-up  or  spin-down,  with  equivalent  energies. 
The  degeneracy  of  the  state  is  therefore  two.  Finally,  if  the  vacancy  is  a  doubly  ionized  acceptor,  the  lowest  energy 
configuration  for  the  two  additional  electrons  is  with  one  spin-up  and  one  spin-down,  with  a  net  state  degeneracy  of 
one.  This  assignment  of  degeneracies  -  one,  two,  and  one  for  the  neutral,  singly  ionized,  and  doubly  ionized  acceptor, 
respectively  -  will  hold  even  if  the  state  is  a  Hubbard  negative-U  state,  although  in  this  case,  the  singly  ionized  state 
will  not  be  occupied. 

A  similar  argtunent  follows  for  the  other  donor  and  acceptor  defect  levels.  In  general,  for  the  II- VI  materials  we 
find  a  degeneracy  of  one  for  the  neutral  defect  state,  two  for  the  singly  ionized  state,  and  one  for  the  doubly  ionized 
state. 


G.  Intrinsic  reaction  constant 

We  are  interested  in  calculating  native  point  defect  densities  at  the  relatively  high  temperatures  at  which  equilibra¬ 
tion  occurs.  It  is  difficult  to  calculate  Kpn  theoretically  because  of  the  difficulty  in  calculating  the  finite-temperature 
band  structure;  in  general  the  band  gap  and  the  conduction  and  valence  band  shapes  ate  all  temperature  dependent. 
An  additional  complexity  in  calculating  Kpn  in  low  x  Hgi.cCdcTe  is  the  nonparabolicity  of  the  conduction  and 
light-hole  bands  near  their  extrema.” 

For  the  purposes  of  evaluating  the  defect  concentrations,  we  have  calculated  Kpn  using  Eq.  (10)  with  ml  =  0.443”; 
an  empirical  relationship  for  the  dependency  of  the  band  gap  on  composition,  x,  and  temperature” 

£,(x,r)  = 

=  -0.302  -I-  1.93z  -  0.810**  -I-  0.832** 

+5.35  X  10-‘‘r(l  -  2*),  (21) 

which  was  fit  for  4.2  <  T  <  300  K;  and  a  linear  dispersion  relationship  for  the  conduction  band  with  a  chosen  to 
yield  good  agreement  with  experimental  values  of  the  intrinsic  carrier  concentrations”  for  T  <  400  K.  We  assume 
that  the  intrinsic  reaction  constant  thus  computed  is  valid  at  temperatures  up  to  655*C,  although  there  have  been  no 
measurements  above  '»400  K  to  substantiate  this  extrapolation.” 


rv.  CALCULATION  OF  DEFECT  FORMATION  FREE  ENERGY:  VIBRATIONAL  CONTRIBUTION 

When  a  defect  is  introduced  into  the  lattice,  the  vibrational  modes  of  the  system  are  modified.  We  must  include  in 
our  calculation  of  the  defect  formation  free  energy  a  term  that  comes  from  modifications  of  the  vibrational  spectrum. 
Most  authors  neglect  this  contribution  to  the  formation  free  energy.  As  we  will  see  in  Section  V,  although  the  electronic 
contribution  to  the  free  energy  is  dominant,  the  vibrational  changes  can  be  significant  and  they  make  a  substantial 
impact  on  the  calculated  magnitude  of  the  defect  concentrations. 

Although  ideally  the  vibrational  contribution  to  the  formation  energy  should  be  calculated  within  LDA  on  the 
same  footing  as  the  static  electronic  contribution,  including  all  of  the  enharmonic  terms,  this  is  a  difficult  and 
computationally  demanding  task.  Instead,  we  take  an  alternative  approach  and  calculate  the  vibrational  spectrum  of 
the  zincblende  lattice  using  Keating’s  valence  force-field  model  for  the  short-range  elastic  interactions.**  Although 
experimental  elastic  constants  were  used,  LDA  theory  actually  predicts  the  elastic  constants  within  10%  for  HgTe 
and  CdTe,  so  we  could  equally  well  have  used  the  calculated  values.  Because  we  are  dealing  with  an  ionic  crystal, 
we  have  included  a  point-charge  model  to  account  for  the  Coulombic  interactions.*^  Unlike  the  valence  force-field 
contributions  to  the  dynamical  matrix,  the  Coulomb  contributions  are  long-range  in  nature  and  induce  a  macroscopic 
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electromagnetic  field,  which  results  in  a  screening  of  the  transverse  optical  phonons.  The  ionic  charge  is  chosen  to 
yield  agreement  with  experiments  for  the  zone  center  splitting  of  the  transverse  and  longitudinal  optical  phonons.  A 
Green’s  function  approach  is  used  to  evaluate  the  lattice-defect  induced  modifications  to  the  phonon  spectrum;  from 
the  perturbed  phonon  density  of  states  the  change  in  the  vibrational  free  energy  can  be  calculated.  Like  the  electronic 
energies,  the  calculations  were  done  for  pure  HgTe  and  are  assumed  applicable  to  defect  calculations  in  Hg0.sCd0.3Te. 
We  expect  that  this  is  a  reasonable  assumption  because  the  elastic  constants  for  HgTe  and  CdTe  are  nearly  identical. 
In  this  paper  we  only  consider  the  vibrational  free  energy  of  the  neutral  defects  and  assume 

FfK  ~  (22) 

for  all  ionization  states.  Details  of  the  calculations  are  given  in  the  Appendix.  Preliminary  results  of  this  model  were 
given  previously.**-^^ 


V.  RESULTS  AND  DISCUSSION 
A.  Formation  free  energies  in  HgTe 

Calculated  defect  formation  electronic  energies  in  HgTe  for  the  defect  reactions  listed  in  Table  I  and  with  Ajt  as 
the  free  mercury  atom,  are  listed  in  Table  III.  For  all  defects,  the  gradient  correction  for  the  16-atom  supercell  is 
calculated  for  the  relaxed  configuration,  as  determined  by  the  LDA  calculation. 

Although  the  total  electronic  formation  energies  listed  in  Table  III  are  important  contributions  to  the  formation 
free  energy,  these  energies  alone  cannot  be  used  to  assess  the  relative  importance  of  the  various  defects  in  the  solid. 
This  is  mostly  due  to  the  free  energy  of  the  atom  in  the  reference  state  (in  excess  of  its  free  atom  electronic  energy) 
that  is  not  included  in  these  electronic  energies,  and  which  is  discussed  in  the  next  section.  This  point  should  be 
obvious  because  we  could  have  just  as  well  defined  our  defect  reactions  with  respect  to  the  tellurium  molecule  in  the 
vapor  phase  and  the  HgTe  solid,  and  obtained  the  corresponding  reaction  energies  that  would  be  quite  different  from 
those  in  Table  III. 

The  calculated  phonon  dispersion  curve  for  HgTe  is  given  in  Fig.  1,  and  is  in  fair  agreement  with  the  experimental 
results.  The  discrepancies  with  the  experimental  curves,  in  particular  near  the  BriUouin  zone  boundary,  can  be 
attributed  for  the  most  part  to  our  neglect  of  long-range  elastic  interactions  in  the  near-neighbor  valence  force-field 
model.^  The  vibrational  entropy  and  energy  contribution  to  the  defect  formation  free  energy  are  calculated  from  the 
density  of  phonon  states;  results  at  500*C  are  given  in  Table  IV.  Eq.  (A20)  can  be  used  to  estimate  the  values  at 
other  (high)  temperatures. 


B.  Defect  reaction  constants  in  Hgo.1Cdo.3Te 


The  reaction  constrmts  for  each  of  the  defect  reactions  listed  in  Table  I  are  calculated  as  a  function  of  temperature, 
where  the  reference  state  was  taken  as  the  monoatomic  mercury  vapor  at  pressure  Pug-  Electronic  energies  for  the 
neutral  defect  formation  energies  are  taken  from  Table  III.  Vibrational  free  energies  are  calculated  using  the  general 
temperature  expression,  as  discussed  in  the  Appendix,  although  for  the  purpose  of  obtaining  an  analytical  expression 
for  the  reaction  constants  with  the  primary  temperature  dependency  explicitly  displayed,  we  have  fit  our  results  for 
500”C  to  the  high-temperature  expression,  Eq.  (A20). 

A  third  contribution  to  the  formation  free  energy  comes  from  the  free  energy  of  the  mercury  in  the  reference  state, 
less  the  electronic  energy  of  the  free  mercury  atom  that  is  contained  in  the  electronic  defect  formation  energies.  For 
the  vapor,  we  must  include  the  translation  free  energy  of  the  atoms.  The  chemical  potential  for  monoatomic  mercury 
in  the  gas  phase  is  given  by  the  standard  expression  from  statistical  mechanics 


HHg  =  kT]n 


(23) 


Combining  these  three  contributions  to  the  reaction  free  energy  we  obtain  the  reaction  constants  for  each  of  the 
defect  reactions.  For  example,  for  the  neutral  mercury  vacancy  in  HgTe,  we  obtain 


“y,.  =  W,!  =  »(i  - 

(1.27  X  10*T-*)exp  . 


(24) 
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where  we  have  taken  g^*  =  1,  -Fvu**  =  -2.09,  0  =  1.48  x  10**  em~*,  and  where  we  have  replaced  txp{-P‘**/kBT) 

by  exp(-(-0.2eV  +  8.9irBT)/ijBr)|6oo«c  =2  1-27  x  10®  T~®.  Reaction  conatants  for  the  eight  native  point  defecU 
considered  in  this  paper  are  listed  in  Table  V. 

The  reaction  constants  for  the  ionized  defects  are  calculated  using  £q.  6  for  acceptors  and  its  generalization  for 
donors.  We  define 


Kx‘'  =  [Xn  =  [X^]^txpixnr  -Ei-  ...El) 

=  Jfx*  }i7exp(zMf  -El-  ...El)  (25) 

and 

Kx>»  =  [•X’**]  =  Kx*  — — *xp(El  +  ...  +  El  -  z^ip).  (26) 

0x* 


C.  Defect  concentrations  In  Hgo.iCdo.jTe 

Gibbs’  phase  rule  tells  us  that  for  a  system  of  three  components  (^4,  B,  and  C)  and  two  phases  (zincblende  solid 
and  vapor)  there  are  three  degrees  of  freedom.  In  evaluating  the  defect  concentrations  in  Hg0.aCd0.3Te  we  have  chosen 
the  temperature,  the  mercury  pressure  and  the  alloy  composition  x  as  these  specified  variables;  the  tellurium 
and  cadmium  pressures,  the  crystal  stoichiometry,  and  the  density  of  the  various  native  point  defects  are  determined 
by  these  conditions. 

The  reaction  constants  in  Table  V  are  evaluated  to  determine  the  concent,  at  ions  of  the  various  native  point  defects 
as  a  frmction  of  temperature  and  pressure.  The  Fermi  energy  is  determined  by  requiring  charge  neutrality.  The 
activation  energies  for  the  shallow  donor  and  acceptor  states  are  taken  to  be  zero;  the  sensitivity  of  our  results  to  this 
assumption  is  discussed  further  below. 

Figs.  2(a)-(c)  show  the  defect  concentrations  at  various  equilibration  temperatures.  Pressure  ranges  are  chosen  so 
as  to  stay  within  the  stability  region  of  the  material.®  At  ^  temperatures  and  pressures  considered,  the  dominant 
defect  is  found  to  be  the  doubly  ionized  mercury  vacancy,  in  agreement  with  previous  interpretations  of  experiments^  ; 
our  result  confirms  the  generally  accepted  experimental  observation  that  the  mercury  vacancy  is  responsible  for  the 
p-type  behavior  of  undoped  HgCdTe  equilibrated  at  high  temperatures 

At  all  temperatures,  the  second  most  dominant  defect  is  found  to  be  the  tellurium  antisite.  The  antisite  concen¬ 
tration  decreases  more  rapidly  with  than  does  the  mercury  vacancy,  and  thus  is  most  important  at  low  mercury 
pressures.  As  does  the  mercury  vacancy,  the  tellurium  antisite  defect  accommodates  excess  tellurium  in  the  lattice, 
and  therefore  its  presence  also  shifts  the  stoichiometry  towards  the  tellurium-rich  side  of  the  phase  diagram. 

The  reason  the  tellurium  antisite  concentrations  is  so  high  deserves  comment.  The  tellurium  antisite  formation 
energy  is  larger  than  that  for  the  mercury  vacancy  by  1  eV  (Table  III).  However,  as  can  be  seen  from  Table  V, 
the  pre-exponential  factor  of  the  reaction  constant  for  the  tellurium  antisite  is  enormous.  The  large  pre-exponential 
factor  results  from  the  large  phase  space  factor  (entropy)  gained  by  creating  two  free  mercury  atoms  compared  to 
that  lost  by  elimination  of  a  formula  unit.  In  contrast,  the  mercury  antisite  density  is  low,  despite  the  fact  that  its 
formation  energy  is  quite  small;  this  again  is  a  result  of  the  pre-exponential  factor  that  in  this  case  is  very  small. 
Thus,  it  is  clear  that  one  must  be  cautious  in  deducing  the  relative  populations  of  the  various  defects  based  on  the 
electronic  contributions  to  the  defect  formation  energies  alone. 

While  the  tellurium  antisite  is  never  found  to  be  the  dominant  defect  controlling  the  doping  under  equilibrium 
conditions,  it  can  introduce  significant  compensation  at  low  mercury  pressures.  Moreover,  the  diffusion  coefficient 
of  the  antisite  ie  expected  to  be  quite  small  because  the  diffusion  of  an  antisite  will  necessarUy  involve  at  least  one 
additional  point  defect,  such  as  the  mercury  vacancy  or  the  tellurium  interstitial.  Thus,  the  tellurium  antisite  may  not 
reach  equilibrium  densities  for  the  times  and  temperatures  corresponding  to  the  Ic  v-temperature  (-^250*0),  high- 
mercury-pressure  anneals  typically  employed  to  reduce  the  mercury  vacancy  density.  If  tellurium  antisite  densities  are 
in  fact  equilibrated  at  a  temperature  at  which  the  antisite  diffusion  effectively  stops  during  cool-down  from  the  growth 
temperature,  then  the  antisites  may  be  frozen-in  at  higher,  nonequilibrium  concentrations.  If  such  a  freezing-in  of 
nonequilibrium  antisites  does  occur,  the  tellurium  antisite  may  in  fact  be  the  "universal”  residual  donor  observed  in 
densities  of  ~  10^®  cm"*  in  materi^  subjected  to  a  low-mercury-pressure,  high-temperature  anneal,  although  it  is 
not  clear  why,  for  example,  the  frozen-in  density  of  antisites  would  be  the  same  for  LPE  material  grown  from  both 
mercury-  and  tellurium-rich  melts.  The  tellurium  antisite  may  also  be  the  defect  responsible  for  the  n-type  carrier 
concentrations  in  as-grown  MBE  material,®  which  is  believed  to  be  grown  on  the  tellurium-rich  side  of  the  phase 
diagram  where  antisite  populations  are  highest. 
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Annealing  strategies  for  reduction  of  the  tellurium  antisite  densities  can  be  developed,  and  may  be  important  if  the 
antisite  is  the  residual  donor.  Consider,  for  example,  a  two-temperature  annealing  process  in  which  a  hrst  anneal  is 
done  under  mercury-saturated  conditions,  but  at  the  lowest  temperature  for  which  the  antisite  is  able  to  equilibrate 
in  reasonable  times.  This  anneal  would  serve  to  lower  the  antisite  densities  as  much  as  possible.  A  second  anneal 
would  be  much  like  that  currently  employed,  -  that  is,  at  ~250*C  under  mercury-saturated  conditions,  -  and  would 
serve  to  anneal  out  the  mercury  vacancies,  leaving  the  antisite  densities  effectively  unchanged. 

At  500*  C  mercury  interstitials  are  present  at  levels  =;  10^°  cm~^  and  at  no  temperatures  are  present  at  levels  high 
enough  to  significantly  compensate  the  mercury  vacancies,  much  less  to  turn  the  material  n-type  under  equilibrium 
conditioiu.  Unlike  the  tellurium  antisites,  the  mercury  interstitials  are  relatively  fast  diffusers,^  and  thus  it  is  unlikely 
that  nonequilibrium  densities  of  interstitials  will  be  frozen-in.  The  mercury  interstitial  densities  that  we  predict  are  in 
quantitative  agreement  with  those  needed  to  explain  the  diffusion  in  Hgi-iCdrTe  in  the  process  simulator  developed 
by  Melendez  and  Helms.^’  There  is  some  uncertainty  in  the  quantitative  predictions  of  the  mercury  interstitial  densities 
reported  here  because  of  the  neglect  of  the  alloy  effects  that  we  expect  to  be  more  significant  than  for  the  mercury 
vacancy  and  tellurium  antisite.  Although  this  correction  will  be  largest  for  the  mercury  interstitial  surrounded  by  four 
mercury  first  neighbors,  it  should  also  be  significant  for  the  interstitial  surrounded  by  four  tellurium  first  neighbors 
because  of  the  six  cation  second-nearest  neighbors,  which  are  only  slightly  more  distant  than  the  first  neighbors. 

The  mercury  antisite  and  the  tellurium  vacancy  and  interstitial  densities  are  all  quite  low,  never  exceeding  10* 
cm~*  at  500*C.  The  corrections  to  these  predicted  densities  may  be  sizable  because  of  alloy  effects,  but  such  corrections 
should  not  significantly  impact  the  densities  of  the  mercury  vacancy  and  tellurium  antisite. 

Fig.  2  shows  the  defect  concentrations  at  the  temperatures  at  which  equilibration  takes  place;  in  Fig.  3  we  show 
the  defect  concentrations  for  material  equilibrated  at  500*C,  then  quench  cooled  to  77*C.  We  have  assumed  that 
the  total  defect  concentrations  are  frozen-in  during  the  quench  -  for  example  [VHg]‘®*“'  =  [V^j]  +  [V^g]  +  [V^g] 
is  constant  -  but  that  the  electrons  and  holes  are  allowed  to  reach  a  new  equilibrium  corresponding  to  the  low 
temperature.  Fig.  4  shows  the  low-temperature  hole  concentrations  for  such  quench-cooled  materials  as  a  function 
of  and  Tanntai,  Compared  with  the  experimental  results  of  Vydyanath.^  The  agreement  of  our  theoretical  results 
with  the  experiments  is  remarkably  good  considering  that  our  calculated  results  are  obtained  almost  entirely  from 
first  principles.  Moreover,  there  is  a  significant  uncertainty  in  the  quenching  efficiency  of  the  experiments;  thus,  exact 
agreement  with  the  experiments  is  not  a  valid  criterion  for  testing  the  accuracy  of  the  theory. 

From  Fig.  4  one  can  see  that  for  the  higher  annealing  temperatures  our  calculations  predict  a  different  dependence 
of  hole  concentration  on  the  mercury  pressure  than  is  experimentally  observed.  The  lower  slopes  for  the  theory  result 
from  our  finding  that  the  material  is  extrinsic  at  the  higher  temperatures.  If  the  material  is  intrinsic  when  equilibrated 
as  the  experiments  indicate,  then  [VHg])it  oc  PjJg ,  as  can  be  obtained  from  Table  V,  with  nr  independent  of  [Viig]. 
However,  if  the  material  is  extrinsic  when  equilibrated  -  that  is,  with  [A*]  =  2[VHg]  -  then  from  the  reaction 

Hg-7,i'g  +  2A*  +  Hg„p,,  (27) 


we  see  that  [V,i'g][A*]*  «  [Vji'g]*  oc  ^Hg. 

The  discrepancies  between  theory  and  experiments  may  be  due  to  a  number  of  factors.  First,  both  uncertainty  in  the 
quenching  efficiencies  and  analysis  of  the  Hall  data  may  account  for  some  of  the  discrepancy.  There  are  also  a  number 
of  uncertainties  in  the  theoretical  calculation  that  may  account  for  the  discrepancies.  These  include  uncertainties  in 
the  electronic  and  vibrational  defect  formation  free  energy,  the  ionization  energies  of  the  defects  (which  were  assumed 
to  be  zero  in  the  above  calculations),  alloy  effects,  and  finally  the  uncertainties  in  the  intrinsic  reaction  constant. 
These  are  discussed  in  turn  below. 

First  it  is  interesting  to  examine  the  sensitivity  of  our  predictions  to  the  accuracy  of  electronic  and  vibrational  defect 
formation  free  energy.  In  Fig.  5  we  have  recalculated  the  77  K  hole  concentrations  as  a  function  of  the  annealing 
temperature  with  the  electronic  contribution  to  the  mercury  vacancy  formation  energy  increased  by  10%,  and  a  rigid 
shift  upwards  in  our  results  by  a  factor  of  2.5;  such  small  changes  result  in  better  agreement  with  experiments.  The 
corrections  to  our  calculated  mercury  vacancy  and  tellurium  antisite  formation  energies  due  to  alloy  effects  have  not 
yet  been  included  in  these  calculations  and  are  expected  to  be  in  the  range  of  tenths  of  an  eV,  as  discussed  in  Section 
lie,  and  thus  may  account  for  the  magnitude  of  correction  used  in  this  example.  Preliminary  estimates  for  the 
corrections  for  going  from  the  16-atom  to  the  32-atom  supercell  are  —O.l  eV  for  the  cation  vacancy  and  the  tellurium 
antisite  as  well.  Finally,  Jahn-Teller  relaxation  energies  have  not  been  included  in  the  present  work  and  they  may 
modify  the  electronic  formation  free  energies.  An  increase  in  the  effective  vibrational  frequencies[(w  in  Eq.  (A20)] 
can  account  for  an  upward  shift  in  the  densities.  Such  an  increase  may  arise  from  differences  between  the  neutral 
and  ionized  defect  vibrational  free  energies  and,  perhaps,  anharmonic  effects  that  may  be  large  at  defects  such  as  the 
vacancy  where  an  atom  is  missing  from  the  lattice.  Thus,  we  see  that  our  calculations  agree  with  the  experimental 
data  approximately  to  within  the  known  uncertainties  of  the  theory. 
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Our  results  are  also  very  sensitive  to  the  intrinsic  reaction  constant,  which  in  turn  depends  sensitively  on  the  band 
structure  and  its  temperature  dependence.  HgCdTe  is  known  to  be  anomalous  in  that  its  band  gap  is  found  to 
increase  with  temperature  at  room  temperature  and  below,  and  although  there  is  no  experimental  information  cm  the 
temperature  dependence  of  the  gap  at  higher  temperatures,  we  have  assumed  that  Eq.  (21)  extrapolates  to  higher 
temperatures.  In  addition,  as  discussed  above,  we  have  assumed  a  parabolic  valence  band,  but  a  linear  variation  of 
the  conduction  band,  with  the  slope  chosen  to  agree  with  the  intrinsic  carrier  concentrations  [with  the  energy  gap 
given  by  Eq.  (21)  at  temperatures  below  400  *C].  While  this  fit  is  quite  good  for  the  temperature  range  over  which 
it  is  fit  (from  77  K  to  400  K),  the  reliability  of  K,n  at  250*  C  and  above  for  which  we  have  presented  our  defect 
density  predictions  is  unknown.  To  demonstrate  the  sensitivity  of  our  results  to  the  intrinsic  reaction  constant  we 
have  calculated  the  77  K  temperature  hole  concentrations  with  the  conduction  band  density  of  states  increased  a 
factor  of  ten  and  have  found,  except  at  the  very  highest  temperatures,  that  the  calculated  hole  concentrations  vary  as 
P^g  indicating  intrinsic  behavior  at  the  annealing  temperature,  and  resulting  in  better  agreement  with  experiments. 
Because  our  results  depend  sensitively  on  the  intrinsic  reaction  constant,  it  is  essential  to  establish  a  reliable  prediction 
of  its  value  at  the  annealing  and  growth  temperatures  where  equilibration  of  the  defect  densities  takes  place.^ 

We  have  suggested  that  nonequilibrium  densities  of  tellurium  antisites  may  be  the  residual  donor,  but  they  may 
also  be  important  SRH  recombination  centers.  It  is  experimentally  observed  that  the  residual  donor  does  not  freeze 
out  even  for  samples  cooled  to  4  K,  and  therefore  its  first  ionization  state  must  resonate  in  the  conduction  band. 
A  SRH  recombination  center  in  Hgi-xCd^Te  with  x  =  0.22  lies  <^25  meV  below  the  conduction  band  edge,  has  a 
larger  capture  cross  section  for  electrons  than  holes,  and  typically  has  a  density  smaller  than,  but  comparable  to,  the 
residual  donor  density.  The  properties  of  the  antisite  are  consistent  with  such  a  level;  it  is  a  donor;  although  we  have 
assumed  here  it  is  a  single  donor,  it  is  likely  that  a  second  donor  level  is  present  in  the  gap  and  may  be  ionized  at 
the  high  processing  temperatures;  and  the  first  ionization  level  may  be  resonant  in  the  conduction  band,  although  we 
are  unable  to  resolve  this  in  our  present  calculations  with  certainty.  A  more  quantitative  prediction  of  the  ionization 
levels  of  the  antisite  is  needed  to  correlate  it  with  a  SRH  center. 

A  technologically  important  step  in  making  IR  detectors  from  Hgi.rCdcTe  is  a  low-temperature  mercury-saturated 
anneal  that  is  done  to  reduce  the  mercury  vacancy  concentrations.  In  Fig.  6  we  show  the  defect  concentrations  for 
material  annealed  at  various  temperatures  along  the  mercury-saturated  side  of  the  phase  diagram.  In  such  mercury- 
saturated  anneals,  if  equilibrium  can  truely  be  reached,  then  the  mercury  vacancies  will  certainly  be  the  dominant 
defect,  with  the  tellurium  antisite  density  being  negligible.  However,  as  discussed  above,  it  is  unlikely  that  equilibrium 
densities  of  antisites  will  be  achieved  at  these  relatively  low  temperatures. 

The  above  analysis  of  defect  concentrations  can  be  repeated  for  a  number  of  diiferent  situations.  For  example,  we  can 
calculate  the  native  defect  densities  with  a  donor  or  acceptor  impurity  present.  At  the  high  growth  temperatures,  the 
impurity  concentrations  would  have  to  be  comparable  to  the  vacancy  concentrations  to  modify  the  high  temperature 
vacancy  concentration.  On  the  other  hand,  even  impurity  densities  of  the  order  of  10**  cm“*  will  alter  the  native 
defect  concentrations  subjected  to  a  low-temparature  post-growth  anneal.  We  can  also  repeat  the  above  calculations 
for  reference  state  other  than  the  mercury  vapor.  For  example,  in  HgCdTe  tellurium  precipitates  are  known  to  form 
as  a  metastable  state  upon  cooling  from  high  growth  temperatures.**  In  the  vicinity  of  a  precipitate,  the  native  defect 
populations  will  be  in  local  equilibrium  with  the  tellurium  solid,  and  defect  concentrations  for  this  reference  state 
can  be  calculated.  Because  this  constitutes  a  nonequilibrium  situation,  on  must  address  diffusion  rates  to  assess  the 
extent  of  the  modified  defect  atmosphere  about  a  precipitate. 


VI.  SUMMARY  AND  CONCLUSIONS 

We  have  made  quantitative  predictions  of  the  native  point  defect  densities  in  Hgo.sCdo.zTe  as  a  function  of  tem¬ 
perature  and  pressure  and  find  good  agreement  with  the  available  experiments.  We  have  substantiated  the  claim 
that  the  primary  defect  is  the  mercury  vacancy,  and  have  identified  the  tellurium  antisite  as  an  important  secondary 
defect.  A  first-principles  approach  was  used  for  most  of  the  quantities  calculated,  with  the  only  significant  empirical 
data  being  those  needed  to  obtain  the  temperature-dependent  intrinsic  reaction  constant. 

Although  we  predict  the  undoped  material  to  be  always  p-type,  refinements  in  our  calculations  may  show  that  the 
antisite  may  dominate  in  the  low-mercury-pressure  region  and  turn  the  material  n-type  by  a  native  defect;  our  current 
accuracy  is  not  sufficient  to  establish  this.  While  most  euineals  of  technological  importance  are  done  under  mercury- 
saturated  conditions  to  reduce  mercury  vacancy  concentrations,  exploration  of  the  tellurium-saturated  region  where 
we  predict  the  tellurium  antisite  densities  become  comparable  to  those  of  those  of  the  mercury  vacancy  may  help 
confirm  the  presence  of  tellurium  antisites. 

A  second  means  to  explore  the  presence  of  tellurium  antisites  and  their  relationship  to  the  residual  donor  is  through  a 
careful  set  of  experiments  using  two  temperature  anneals,  as  discussed  above.  Because  we  do  not  know  the  temperature 
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at  which  diffusion  of  tellurium  effectively  stops,  the  temperature  of  the  first  anneal  would  have  to  be  varied,  as  would 
the  annealing  time;  the  mercury  pressure  could  also  be  varied,  although  mercury-saturated  conditions  are  those  one 
would  eventually  want  to  employ.  The  identification  of  the  tellurium  antisite  as  the  residual  donor  can  be  made  if  the 
donor  densities  in  identically  grown  material  were  found  to  differ  after  the  second  anneal  (using  the  standard  conditions 
for  a  mercury-saturated  low-temperature  anneal)  depending  on  the  conditions  of  the  first  anneal.  A  quantitative 
analysis  of  this  experiment  would  be  quite  difficult  because,  in  addition  to  uncertainties  in  the  temperature  at  which 
the  teUurium  antisite  equilibration  stops,  if  the  tellurium  antisite  diffuses  via  a  vacancy  mechanism,  the  diffusion  of 
the  antisite  will  depend  on  the  concentration  of  mercury  vacancies  present  during  the  first  anneal. 

It  would  also  be  useful  to  perform  high-temperature  annealing  as  was  done  by  Vydyanath,^  but  instead  of  quenching 
to  77  K  and  having  to  address  the  issue  of  quenching  efficiency,  follow  the  anneals  by  Hall  analysis  at  the  anneal 
temperature.  Recently,  an  attempt  at  such  an  experiment  was  made  by  Wienecke  et  a/.,^‘  although  analysis  of  such 
an  experiment  requires  knowledge  of  the  high-temperature  intrinsic  reaction  constant. 
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APPENDIX:  VIBRATIONAL  ENTROPY  OF  POINT  DEFECTS  IN  SEMICONDUCTORS 
As  discussed  by  Keating,^*  we  assume  the  elastic  energy  for  a  zincblende  structure  can  be  written  as 


where  i  and  i'  sum  over  all  bonds,  A(fifii)  =  rf<,  with  f  and  f°  the  bond  vectors  connecting  adjacent  atoms 

in  the  distorted  and  equilibrium  lattices,  respectively.  For  first-neighbor  interactions,  the  sum  in  the  second  term 
runs  over  only  those  bonds  connected  to  a  common  atom. 

For  a  nonionic  material,  the  normal  modes  for  the  displacement  of  the  atoms  are  determined  by  solving  the  equation 
of  motion  for  the  lattice  cast  in  the  usual  manner  in  terms  of  the  three-dimensional  eigenvalue  equation 


w’tI  =  D*(k)u. 


Here,  D*(k)  is  the  (elastic)  dynamical  matrix  and 


“■'(sO 


is  the  polarization  vector  of  the  normal  modes,  where  u,-  is  the  displacement  vector  of  the  atom.  For  the  zincblende 
lattice  there  are  two  atoms  per  unit  ceil,  so  i=l  or  2. 

In  terms  of  the  elastic  constants,  Cn  and  Cu,  the  dynamical  matrix  for  the  valence  force  field  model  is  given  by 


J 

+  v/3mVm, 


;&s(i:)  +  - 

adc 


7d{Cii-Ci,)i(k) 


1  T 

Tsmj 


where  d  is  the  equilibrium  bond  length,  mi  and  mj  are  the  masses  of  the  two  atoms  in  the  unit  cell,  I  is  the  3x3 
unit  matrix,  and  S  is  given  by 


/  *1  *4  *3 

=  I  54  5i  52 
\  53  52  5i 
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*i(f)  =  -»■  e‘"*  e'*^*  + 

(A«a) 

s,(f)  =  e*"*  +  e*"*  -  e‘”»  - 

(A6b) 

j8(jP)  =  e‘*^‘  -  e‘**  + 

(A6c) 

S4(£)  =  e**^‘  — 

(A6d) 

The  ue  the  vectors  connecting  atom  1  to  atom  2  in  the  unit  cell  and  are  given  by  di  s  dj  =  }[lil], 

da  s  end  ^4  =  f  [III],  where  a  is  the  lattice  constant. 

While  short-range  elastic  forces  are  described  within  the  valence  force  field  model,  in  crystal  with  an  ionic  contri¬ 
bution  to  the  bonding.  Coulomb  interaction  must  also  be  included  in  the  dynamical  matrix.  The  long-range  nature  of 
the  Coulomb  interaction  complicates  the  problem  considerably.  Using  a  pair-wise  point  charge  model  of  the  Coulomb 
interaction,  the  Coulomb  dynamical  matrix  is  given  by^^: 


(A7) 


where 


X  E,<  ffac  (VP(i(l,K)  -  *(/',«'))) 

exp  (=1^) 
exp(i(5(f(«)-fK))). 


(A8) 


In  the  above  equation  I  and  k  label  the  unit  cell  and  basis  atoms,  a  refers  to  the  cartesian  component,  g,  is  the  effective 
charge,  S(k,  1)  s  f  (k)  +  £{l)  is  the  postition  vector  of  the  «c'*  atom  in  the  unit  cell,  G  are  the  reciprocal  lattice 
vectors,  vo  is  the  unit  cell  volume,  cq  is  the  permitivity  constant,  and  P  is  a  (numerically  determined)  measure  of  the 
Gaussian  charge  distribution  used  in  the  Ewald  summation.  Ha^„>{y)  is  given  by  the  integro-differential  expression 


The  full  dynamical  matrix  is  given  by  the  sum  of  the  Coulomb  and  elastic  contributions 


D(*)  =  D'(*)-J-D*(*). 


(A9) 


(AlO) 


In  the  present  problem  we  are  interested  in  calculating  the  change  in  the  vibrational  free  energy  of  the  crystal  lattice 
due  to  the  creation  of  a  defect.  This  is  done  using  the  Green’s  function,  which  is  most  conveniently  calculated  in  terms 
of  the  density  of  states  of  the  phonon  system.  For  the  ideal  crystal  without  a  defect,  a  Brillouin  zone  integration  is 
done  to  calculate  the  phonon  density-of-states  matrix,  0,  from  which  the  Green’s  function  can  be  calculated  via 


gV)= 


(All) 


where  the  singularity  in  the  integral  has  been  explicitly  removed.  The  change  in  the  total  density  of  states  when  a 
defect  is  introduced  into  the  crystal  can  be  deduced  from  Dyson’s  equation  to  obtain 


A0(u;*)  =  ^Im^ln  (det(l  -  G°(w2)V)) 


(A12) 
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where  V  ia  the  perturbation  potential.  In  the  present  case  we  use  a  strictly  site-diagonal  perturbation  potential 
corresponding  to  the  mass  change  due  to  the  introduction  of  an  isdated  defect. 

The  partition  function  for  the  phonon  system  in  the  sincblende  lattice  is  given  by 


(A13) 

where  Ui  are  the  normal  modes  of  the  system.  The  total  vibrational  entropy  of  the  system  is  obtained  from  s 

^{kaTlaZ),  which  gives 

(A14) 

We  convert  the  sum  to  an  integral  by  the  replacement 

5^  — *  J  p(u)du  =  J  2p(u^)udu 

(A15) 

where  p(w^)  =  TrO(u^).  We  are  interested  in  the  change  in  the  vibrational  entropy  upon  formation  of  a  defect,  in 
which  case  we  replace  p  by  Ap  to  obtain 

S'"  =  St,  jf"  M-’)  (*) 

—In  ^2sinh  (s^)))"'^-  (*“> 

Similarly  for  the  vibrational  energy,  we  use  the  relationship  IP**  s  -g^f;y^ln2  to  obtain 

cr»  =  s_("Mw’)  (1.0.1.  (AIT) 

with  the  change  in  the  vibrational  free  energy  F***  =  IP**  —  TS***. 

While  the  calculation  of  the  vibrational  terms  in  this  paper  was  done  using  the  general  expressions  above,  it  is 
interesting  to  examine  the  expression  for  the  free  energy  in  the  high  temperature  limit,  that  is  appropriate  for  high 
growth  and  processing  temperatures,  and  to  examine  the  explicit  temperature  dependence  of  this  term.  In  the  high 
temperature  limit  hwmac  <  ^bT  uid  P'**  reduces  to 

«  2ibflr  In  p{u)u(L).  (A18) 


This  integral  can  be  shown  to  be  equal  to 

(A19) 

where  is  the  number  of  phonon  modes  created  oi  -li^troyed  in  the  defect  reaction  of  interest  ,  and  w  is  an  appropri¬ 
ately  weighted  frequency.  What  enters  the  calculation  of  the  defect  concentrations  is  txp{—P**/kBT)  which  reduces 
to 


(A20) 


in  the  high  temperature  limit.  For  the  defect  reactions  in  Table  I  we  obtain 
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-3  Vi4 
3  Vt. 

nm=:6  HgT, 

-«  Ten, 

+3  Hg/ 

-3  Te/ 

(A21) 

We  will  UM  this  eimple  power-law  dependence  of  Eq.  (A20)  to  extract  a  tinaiple  power-law  temperature  dependence 
of  the  reaction  conatanta. 
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FIG.  1.  Calculated  phonon  diaperaion  curve  for  HgTe 


FIG.  2.  Defect  concentration!  at  annealing  temperature  for  material  annealed  at  (a)  250*C,  (b)  S00*C,  and  (c)  65S*C. 


FIG.  3.  Defect  concentration!  for  500*C  anneal  after  quench  cooling  to  77  K. 

FIG.  4.  Hole  concentration!  at  77  K  for  material  equilibrated  at  varioua  high-temperature  annealing  conditiona.  Experi¬ 
mental  reanlta  taken  from  Ref.  (1). 

FIG.  5.  Hole  concentration!  at  77  K  for  material  equilibrated  at  varioua  high-temperature  annealing  condition!,  calculated 
with  the  electronic  formation  energy  for  the  mercury  vacancy  increaaed  by  10%  and  a  rigid  upward  ahift  of  the  hole  concentration! 
by  a  factor  of  2.2.  Thia  figure  ia  meant  to  demonatrate  the  aenaitivity  of  our  reaulta  to  amall  changea  in  our  calculated  parameter!. 
Aa  diacuaaed  in  the  text,  aimilar  qualitative  changea  in  our  reaulta  can  be  aeen  by  modifying  the  intrinaic  reaction  conatant. 

FIG.  6.  Defect  concentration!  for  material  annealed  along  the  mercury-rich  aide  of  the  atability  region  at  (a)  the  annealing 
temperature  and  (b)  after  quenching  to  77  K. 

TABLE  I.  Defect  reactione  conaidered  for  compound  AC.  Notation  ia  aa  foUowa;  The  primary  aymbol  refera  to  the  apedea, 
the  aubacript  refera  to  the  aite  that  the  apeciea  occupiea,  with  no  aubacript  indicating  that  the  apedea  ia  occupying  ita  uaual 
lattice  aite.  V  correaponda  to  a  vacancy,  I  an  interatitial,  R  aome  external  reference  atate,  and  X  a  generic  defect.  Two  typea 
of  interatitiala  are  conaidered,  both  occupying  tetrahedral  ntea,  the  firat  aurrounded  by  four  cation  neaieat  neighbor!,  7c,  the 
aecond  aurrounded  by  four  anion  neareat  neighbor!,  Ia-  Following  the  notation  of  Kroger,**  an  x  anperacript  correaponda  to  a 
neutral  apedea,  a  prime  to  a  negatively  charged  apedea,  a  bullet  to  a  poaitively  charged  apeciea,  and  e'  and  h*  are  an  electron 
and  hole,  reapectively. 
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TABLE  n.  Bnlk  cohanve  tncigiea  with  ud  without  the  giadieat  conectioiu  (GC)  to  the  locel  deueity  (LD),  end  companeou 


with  experiment. 


Compound 

LD 

Coheaive  energy  (eV/bond) 

GC 

Experiment 

$i 

2.58 

2.31 

2.32 

Ge 

2.22 

1.88 

1.94 

AlP 

2.35 

2.05 

2.13 

AlAi 

2.17 

1.85 

1.69 

AlSb 

1.91 

1.61 

1.76 

GuP 

2.08 

1.76 

1.78 

GuAi 

1.91 

1.56 

1.63 

1  GuSb 

1.70 

1.36 

1.48 

InP 

1.89 

1.56 

1.74 

In  As 

1.77 

1.42 

1.55 

InSb 

1.80 

1.26 

1.40 

ZnS 

1.82 

1.53 

1.59 

ZnSe 

1.64 

1.35 

1.29 

ZnTe 

1.43 

1.15 

1.20 

CdTe 

1.33 

1.04 

1.10 

HgS 

1.29 

0.94 

1.02 

HgSe 

1.19 

0.84 

0.85 

BgTe 

1.09 

0.76 

0.81 

TABLE  III.  Neutrul  native  defect  formation  energies  for  HgTe  corresponding  to  reactions  in  Thble  I,  where  A  is  mercury,  B 
is  tellurium,  and  the  free  atom  is  used  as  the  reference  state  Ajt.  Local-density  (LD)  calculations  were  done  using  a  16-atom 
supercell,  unless  otherwise  noted.  Gradient  correction  (GC)  energies  are  discusiKd  in  the  text.  The  most  important  ionisation 
states  of  the  native  point  defects  are  also  given. 

Defect 

Energy  (eV) 

Ionisation 

State 

LD 

Rel&x&tion* 

cc 

Total 

Vh, 

2.83 

-0.05 

-0.69 

2.09 

douUe  acceptor* 

0.99 

-0.01 

0.49 

1.47 

shallow  donor 

Hg-p, 

-1.00 

0 

1.46 

0.46 

deep  acceptor 

TeHj 

4.85 

-0.19 

-1.54 

3.12 

shallow  donor 

0.75 

-0.24 

0.70 

1.21 

shallow  donor 

HglT. 

0.81 

-0.31 

0.62 

1.12 

shallow  donor 

Te/K, 

4.78 

-0.57 

-0.83 

3.38 

shallow  donor 

Te/r. 

5.17 

-0.84 

-0.96 

3.37 

shallow  donor 

‘Aasumed  to  be  u  negutive-U  center 


TABLE  V.  Reaction  conatante  for  nential  defecta  in  Ego.tCdo.aTe  correaponding  to  the  defect  leactiona  in  Table  I.  The 
vibrational  contribution  ia  calculated  at  S00*C,  and  fit  to  the  high-temperatnre  poi*er-law  dependence  [Eq.  A19]  eo  aa  to  ahow 
the  explicit  temperature  dependence.  All  calculationa  in  the  paper  were  done  uaing  reaction  conatanta  with  the  more  exact 
expreaaion  for  the  vibrational  free  energiea. 
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Dr.  Marcy  Berding 
Physical  Electronics  Lab. 

SRI  International,  410-33 
Menlo  Park,  CA  94025 

Febuary  2,  1994 

Dr.  Anis  Husain 
ARPA 

3701  N.  Fedrfax  Drive 
Arlington,  VA  22203 

Dear  Dr.  Husain 

In  response  to  your  suggestions  at  the  meeting  with  Arden  Sher  last  month,  we  have 
spoken  to  both  Dr.  Diego  Olego  at  Philips  Lab  and  Dr.  Michael  Haase  at  3M  about 
their  work  in  ZnSe-based  laser  technology.  These  conversations  were  extremely  useful 
and  both  groups  were  interested  in  having  us  pursue  solutions  to  a  number  of  their 
problems. 

As  a  result  of  these  conversations,  we  were  able  to  identify  several  areas  in  which  our 
capabilities  could  be  applied  to  current  problems  in  ZnSe-based  laser  technology.  These 
areas  are  summarized  below: 

DOPING 

Doping  is  still  an  area  of  great  interest  and  remains  a  performance  limiter.  Although 
some  success  has  been  achieved  with  nitrogen,  what  is  limiting  the  whievable  p-type 
doping  levels  is  still  a  mystery  that  both  groups  are  anxious  to  have  solved.  Some  curi¬ 
ous  effects  have  also  been  seen,  such  as  a  reduction  in  the  p-type  conductivity  directly 
under  ZnTe  contacts.  While  both  groups  were  less  optimistic  about  the  prospects  of 
achieving  high  levels  of  p-doping  with  phosphorus  or  arsenic,  they  both  thought  it 
should  be  an  area  of  continued  research. 

We  are  currently  working  on  the  issue  of  doping  under  om  ARPA  sponsorship.  The 
problem  of  native  defects  is  central  to  answering  the  compensation  question.  We  will 
also  be  looking  at  solubility  limits  of  various  dopants.  Our  proposed  work  on  dopants 
was  of  interest  to  both  experimental  groups,  because  it  will  provide  guidance  about 
which  dopant  is  most  likely  to  succeed,  and  an  opportunity  for  cost  reduction.  Un¬ 
fortunately,  ui  our  current  rate  of  funding,  it  will  be  over  a  year  imtil  we  finish  these 
calculations. 
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The  efficiency  of  dopants  is  often  related  to  complex  formation.  Both  groups  find  there 
is  more  nitrogen  incorporation  than  there  are  active  acceptors.  In  other  materials  this 
is  often  the  result  of  complex  formation,  either  between  impurity  atoms  or  between  an 
impurity  and  a  native  point  defect.  We  are  currently  calculating  complex  formation 
energies  in  HgCdTe  and  heavily  arsenic-doped  silicon.  Our  current  program  on  ZnSe 
does  not  include  work  on  complexes. 

There  a  number  of  extensions  we  csm  make  to  our  work  to  study  the  doping  near 
interfaces  such  as  ZnSe/ZnTe;  for  example,  once  we  have  completed  our  current  work, 
it  would  be  straightforward  to  extend  our  calculations  to  include  band  bending  at  the 
interface,  and  to  predict  the  resulting  modifications  in  the  equilibrium  doping  eflSciency. 
We  can  also  explore  nonequilibrium  impurity  and  defect  densities  that  may  result  when 
growth  occurs  at  temperatures  below  which  diffusion  of  the  relevant  species  stops.  We 
also  have  tools  in  hand  to  consider  modifications  to  the  bulk  doping  that  results  from 
alloying  ZnSe  with  other  elements  such  as  sxilfur,  cadmium,  and  tellurium. 

DARK  LINE  DEFECTS 

This  was  the  area  of  most  interest  to  both  groups,  since  dark  line  defects  are  limiting 
device  lifetimes,  preventing  products  from  being  brought  to  market.  While  the  dark 
line  defects  appear  to  be  related  to  dislocations  and  perhaps  native  point  defects,  the 
full  relationship  is  still  not  understood.  As  devices  degrade,  the  dark  areas  grow, 
but  it  is  not  obvious  whether  it  is  dislocation  multiplication  or  modifications  to  the 
defect  atmosphere  around  existing  dislocations  that  result  in  the  device  failure.  There 
also  seems  to  be  a  correlation  of  device  failure  with  stacking  faults  originating  at  the 
ZnSe/GaAs  interface,  and  removal  of  these  faults  is  of  critical  importance.  While  the 
propagation  of  dislocations  through  the  active  QW  structure  is  correlated  with  failure 
of  the  devices,  the  means  by  which  the  devices  fail  is  still  not  clear;  that  is,  do  the 
dislocations  short  the  jvmction,  or  do  they  just  act  as  recombination  sites. 

The  area  of  the  impact  of  dislocations  on  device  performance  is  an  area  we  have  been 
just  recently  working  in  the  context  of  HgCdTe  double-layer  heterojunction  devices. 
There  we  examined  the  impact  of  both  the  dislocation  strain  field  (and  the  resulting 
pieozoelectric  charge)  and  core  charges.  This  theory  can  be  easily  extended  to  ZnSe- 
based  laser  structures.  In  addition,  this  work  can  be  extended  to  examine  the  nature  of 
the  material  in  the  vicinity  of  the  dislocation,  for  example  native  defect  and  impurity 
populations,  and  their  correlation  with  the  observed  growth  in  the  iffea  of  dark  line 
defects.  While  the  properties  of  stacking  faults  -  for  example,  their  energy  -  can  be 
examined  using  our  tools,  most  of  the  technological  concern  is  in  elimination  of  the 
stacking  faults,  which  is  closely  coupled  to  the  MBE  growth  process,  discussed  further 
below. 
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BAND  STRUCTURE  PROPERTIES 

While  many  ZnSe-based  alloys  are  used  in  current  ZnSe-based  lasers,  very  little  is 
known  about  the  band  offsets  among  the  constituents.  We  were  told  that  knowledge  of 
these  offsets  would  greatly  help  in  device  analysis  and  design.  In  addition,  more  infor¬ 
mation  on  the  band  gap  variation  with  composition  in  alloys  such  as  ZnSSe,  ZnCdSe, 
ZnSeTe,  and  ZnMgSSE  would  be  extremely  useful.  The  area  of  valence-band  offsets  is 
also  of  interest  for  the  new  alloy  system,  ZnSeTe,  being  explored  for  its  ease  of  p-type 
doping  and  ability  to  make  ohmic  contacts 

Once  again,  this  is  an  area  in  which  our  alloy  theory  group  has  had  extensive  experi¬ 
ence.  We  have  calciilated  the  band  structure  of  most  of  the  group  IV,  III-V,  and  II- VI 
compovmds  and  many  of  their  alloys.  Otir  results  have  been  in  quantitative  agreement 
with  experiment.  We  have  also  applied  these  techniques  to  new  alloy  systems,  such  as 
HgZnTe,  and  our  results  were  later  substantiated  by  experiments.  Accurate  calciila- 
tions  of  band  offsets  among  the  II-VI  compounds  and  their  alloys  are  also  possible  using 
a  combination  of  our  techniques:  the  first-principles  full-potential  linearized  muffin-tin 
orbital  and  the  tight-binding-based  coherent  potential  approximation  (CPA)  method. 

MOLECULAR  BEAM  EPITAXY  (MBE) 

The  major  successes  to  date  in  the  ZnSe-based  laser  technology  have  been  with  mate¬ 
rials  grown  by  MBE  and  doped  with  nitrogen.  While  great  progress  has  been  made,  a 
nvimber  of  issues  relating  to  growth  must  still  be  addressed.  One  of  these  relates  to  the 
failure  mechanism  and  its  relationship  to  stacking  faults  at  the  III-V/  II-VI  interface. 
While  MBE  growth  parameters  have  been  optimized  to  eliminate  the  formation  of  in¬ 
terfacial  phases  such  as  Ga2Se3,  elimination  of  stacking  faults  still  remains  an  issue. 
Issues  of  doping  also  closely  relate  to  the  growth  process,  as  is  well  known  from  the 
need  to  use  a  plasma  source  for  nitrogen  cracking  to  obtain  the  current  nitrogen  doping 
levels. 

The  current  limitation  in  the  p-doping  levels  achievable  may  be  related  to  the  growth 
itself,  rather  than  a  fundamental  solubility  limitation  or  native  defect.  In  our  past  work 
we  have  demonstrated  that  in  some  wide-gap  ionic  materieds,  atoms  on  the  surface 
have  a  mutually  repulsive  interaction,  resiilting  in  ordered  arrays  of  surface  atoms 
and  vacancies.  The  occurrence  of  this  ordering  depends  on  the  ciirface  orientation, 
substrate  temperature,  effective  atomic  fluxes  and  surface  stabilization,  and  surface 
coverage.  Substantial  ordering  may  occur  for  only  the  cation  or  anion  layer.  When 
present,  this  ordering  will  modify  the  nature  of  nonisoelectronic  impurity  and  alloy 
constituent  incorporation  into  the  layer.  This  theory  predicts  the  observed  ordering  in 
MBE-grown  GaAlAs,  and  is  currently  being  experimentally  tested  for  silicon  impurity 
incorporation  in  GaAs. 
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PROPOSED  TASK  LIST 

Based  on  our  discussions  with  the  experimental  groups,  very  briefly  outlined  above,  we 
have  developed  a  task  list  for  an  extension  of  our  current  program  on  ZnSe.  Follovring 
the  task  list  is  an  estimate  of  the  man-years  required  for  its  accomplishment.  We  have 
not  provided  a  description  of  the  methods  we  will  employ  in  carrying  out  these  tasks, 
nor  references  to  our  prior  publications  in  the  relevant  fields.  However,  in  all  of  the 
proposed  tasks  we  have  from  moderate  to  extensive  relevant  previous  work. 

TASK  1 

Devise  means  to  improve  p-doping  in  ZnSe-based  materials,  increasing  dopant  incor¬ 
poration  and  efiiciency. 

1.1  Extend  our  current  program  to  calculate  the  equilibrium  solubility  of  isolated  im¬ 
purities  for  the  most  promising  dopants  in  ZnSe,  including  nitrogen,  phosphorus  and 
arsenic.  Determine  if  currently  realized  limits  for  active  nitrogen  incorporation  are  due 
to  equilibrium  solubility  effects  and,  if  not,  what  are  those  limits.  Determine  if  any 
impurity  looks  more  promising  than  nitrogen. 

1.2  Calculate  the  density  of  defect  pairs,  including  impurity-impurity  and  impvirity- 
defect  pairs,  with  the  aim  of  determining  if  pairing  is  limiting  the  doping  efficiency. 
The  impurities  to  be  considered  depend  on  the  predicted  solubilities  in  Task  1.1. 

1.3  Based  on  results  from  Tasks  1.1  and  1.2,  devise  a  strategy  to  improve  the  doping 
levels  obtainable  in  ZnSe-based  materials. 

1.4  Examine  the  nature  of  the  ZnSe/ZnTe  interface  to  seek  an  explanation  for  the 
reduced  solubility  of  nitrogen  there.  Propose  means  by  which  to  increase  the  local 
solubility. 

1.5  Consider  modifications  to  all  of  the  above  based  on  the  fact  that  the  material  may 
be  in  quasi-equilibrium  during  and  following  the  low-temperature  MBE  growth. 

TASK  2 

Develop  means  by  which  to  increase  the  operational  life  of  ZnSe-based  LEDs  and 
Lasers. 

2.1  Extend  our  theory  of  dislocations  and  their  impact  on  HgCdTe  devices,  to  ZnSe- 
based  laser  structures.  Determine  the  means  by  which  dislocations  impact  device  op¬ 
eration  and  correlate  dislocation  density  with  device  efficiency. 

2.2  Determine  impiirity  and  native  defect  atmospheres  about  dislocations  due  to  both 
Coulomb  and  strain  fields.  Consider  how  these  atmospheres  change  during  the  op¬ 
erational  lifetime  of  devices,  including  effects  of  nonradiative-recombination-induced 
bond-breaking  near  the  p-n  junctions  and  thermal  heating. 
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2.3  Based  on  the  results  of  Task  4,  examine  mechanisms  responsible  for  causing  stacking 
faults  in  MBE  growth  at  the  ZnSe/GaAs  interface,  and  devise  strategies  to  eliminate 
them. 

TASKS 

Provide  parameters  needed  for  design  of  present  and  future  ZnSe-based  Lasers. 

3.1  Calcxilate  the  band  structures  of  ZnSe  and  related  alloys.  Calculate  the  band-gap 
bowing  as  a  function  of  alloy  composition.  These  calculations  will  also  provide  the 
wave  functions  needed  for  other  tasks  listed  below. 

3.2  Predict  the  band  offsets  among  the  alloys  and  compounds  relevant  to  the  ZnSe- 
based  laser  technology,  as  specified  by  the  experimental  groups. 

3.3  Calculate  the  nonradiative  lifetimes  for  the  relevant  mechanisms  of  Auger,  and 
Shockley- Read  for  impurities  and  native  point  defects  present  in  significant  concentra¬ 
tions. 

TASK  4 

Optimize  MBE  growth  parameters  for  incorporation  of  active  dopants  and  for  mini¬ 
mization  of  structtiral  defects. 

4.1  Apply  our  ctirrent  models  of  surfaces  to  ZnSe  and  related  materials  to  predict  the 
nature  of  the  surfaces  during  MBE  growth.  Determine  the  natme  of  the  impurity-host 
interactions  on  the  surface  and  predict  optimal  parameters  for  active  dopant  incorpo¬ 
ration. 

4.2  Develop  a  full  MBE  growth  simxilator  based  on  our  current  first-principles  molec¬ 
ular  dynamics  methods.  Apply  the  simulator  to  the  problem  of  the  GaAs/ZnSe-based 
material  interface  to  predict  the  occinrence  of  stacking  faults;  optimize  growth  pa¬ 
rameters  to  minimize  stacking  fault  densities.  Exercise  the  model  in  numerous  other 
applications  relevant  to  the  growth  of  ZnSe-based  laser  devices. 

Assuming  a  3-year  effort,  we  have  estimated  1.25  man-years  to  accomplish  Twk  1, 
0.75  man-year  for  Task  2,  0.5  man-year  for  Task  3,  and  1.5  man-years  for  Task  4. 
While  most  of  the  work  is  devoted  to  the  ZnSe-based  systems,  the  development  of  an 
MBE  simulator  under  Task  4  would  be  applicable  to  many  other  systems  in  addition 
to  ZnSe. 

Tasks  1  and  4  axe  especially  computationally  intensive,  and  will  require  the  purchase 
of  a  workstation  with  about  128  MBytes  of  memory.  We  estimate  this  would  cost 
approximately  $50,000. 

Although  neither  Dr.  Haase  nor  Dr.  Olego  informed  us  that  our  conversations  were 
confidentied,  we  hesitate  to  send  them  this  letter  because  it  contains  information  ex¬ 
tracted  firom  both  conversations.  However,  we  certsdnly  will  be  comfortable  if  you 
choose  to  send  copies  of  this  letter  along  to  them  for  feedback. 
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We  are  very  excited  about  getting  a  start  on  the  tasks  outlined  above.  We  feel  these 
tasks  are  in  areas  in  which  we  can  effectively  collaborate  with  the  experiment  teams 
and  make  a  substantial  impact  on  the  technology.  We  will  call  you  in  few  weeks  to 
obtain  your  feedback  on  this  proposed  program. 

Sincerely, 

Marcy  Herding 
Senior  Research  Physicist 


Arden  Sher 
Associate  Director, 

Physical  Electronics  Laboratory 


MB/mb 


P.  S.  We  have  just  receive  a  letter  from  Dr.  Atsuko  Ebina  of  Tohoku  University  who 
is  working  on  ZnSeSe  and  ZnSeTe  and  who  has  just  obtained  some  very  interesting 
results  on  the  band  structure  of  these  alloys  as  a  function  of  alloy  composition.  A 
copy  of  her  letter  and  a  couple  of  the  numerous  figures  she  sent  us  are  attached.  She 
will  be  just  down  the  road  at  Stanford  University  this  summer  working  with  Professor 
Spicer  and  would  like  to  collaborate  with  xis.  Because  we  have  done  extensive  work 
in  the  theory  of  alloy  ordering  and  phase  transitions  in  the  semiconductor  alloys,  we 
know  that  we  could  have  an  effective  and  fruitful  collaboration  with  her.  The  resiilts  of 
such  an  effort  would  greatly  benefit  the  U.S.  experimental  groups.  As  such,  we  would 
like  to  include  an  additional  task  to  collaborate  with  Dr.  Ebina  as  well  as  the  U.S. 
experimental  groups  and  to  apply  oiu:  alloy  theory  of  the  phase  transition  to  the  ZnSSe 
and  ZnSeTe  systems.  We  expect  this  would  require  of  the  order  of  0.5  manyear. 
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